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Abstract 


The combined capacity expansion, production- 
inventory problem is modeled as a linear, integer 
program and then extended to also allow capacity de- 
crease. The models assume constant returns to scale 
in the production function of a firm which must meet 
deterministic demands for a single product over N 
periods, with no backordering, at minimum cost. 

A linear transformation is used to obtain 
equivalent forms of the models which are then decompos- 
ed into fixed cost and variable cost parts.  Enumer- 
ation is done on the fixed cost variables and transpor- 
tation sub-problems are solved with the remaining var- 
fables. Special demand and cost structures are inves- 
tigated, one of which leads to formulation of the en- 
tire problem as a network. 

Extensions are made to include piecewise- 
linear, convex production costs and lumpiness in ca- 
pacity acquisition. "A reinterpretation of Che models 
to handle the employment smoothing problem is also pre- 
sented. Finally, computational results of a Fortran IV 


coding of the basic algorithms are presented. 
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Chapter I 


I.1 Introduction 


The quantitative approach to the solution 
of problems in the realm of managerial economics can 
be traced back to the efforts of the early industrial 
engineers in the United States at the end of the nine- 
teenth century. These men, most notably Frederick W. 
Taylor, pioneered the industrial engineering concepts 
such as Job description, time-motion studies, etc., 
which had such a great influence on the development 
of the assembly line in industry during the first 
half of the twentieth century. 

with World War II, however, an entirely new 
set of tools was refined and applied by the early prac- 
titioners in the infant field of Operations Research. 
Though most of the techniques had their origin in 
military applications, they gradually found their way 
into the industrial sphere in the peace-time economic 
expansion resulting from the war. In the intervening 
decades, the complexity of the practical problems 
Solved has moved apace with the rapid advances in theory 
typical of a new field. The logical result is that an 


increasing amount of effort has been put into inte- 





gration of past techniques using newer theory in an 
attempt to provide more realistic approaches to problems 
with broader horizons. 

The present thesis makes such an attempt by 
offering models of the combined production-inventory 
and capacity expansion problem and algorithms which 
give all integer solutions to these models. Dis- 
cussion of the importance of the combined problem and 
the integral approach are given later in this chapter. 
This is the first integer approach to the combined 
problem known to exist. 

The basic technique used will be branch and 
bound with network sub-problems. A unique feature 
of the algorithm ís the capability of starting it 
at any point in the enumeration. This is especially 
meaningful when linear programming is used to give a 
universal, and often non-integer, optimal solution. 
The enumeration can be started from a suitable point 
with respect to the linear programming optimum to 
determine the existence of any optimal, rounded-off 
solutions rather quickly. 

The famous economic problem of "lumpiness" 
in capacity acquisition as well as multi-shift, and 
overtime, operation of the production facility will 


be treated. A discussion of the relationship between 





the traditional capacity expansion and employment 
smoothing literatures and re-interpretation of the re- 
sults in terms of employment smoothing will also be 


given. 


I.2 Problem Description 


A useful paradigm is a company with a single 
producing facility and a homogeneous product faced 
with a schedule of deterministic demands over a set of 
discrete periods to some pre-set horizon, These de- 
mands must be met on time from production or existing 
inventory at minimum possible cost. 

The capacity of the production facility is 
defined as its output in units per time period when 
equípment is being run at its maximum efficient rate 
with a single shift not working overtime. The models 
will be extended to allow multiple shifts and over- 
time work, but note that these additions are not ger- 
mane to the basic definition of capacity. 

The company's basic problem, assuming that 
cumulative demand never exceeds cumulative capacity 
(where cumulative capacity now includes all regular 
and overtime production), is to find a production- 


inventory schedule to meet demand at minimum 
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4, 
cost.1 References to this problem are ubiquitious in 
both the operations research literature and the ec- 
onomics literature, and will not be discussed in this 
work. 

If, however, cumulative demand does exceed 
cumulative capacity within the planning horizon, then 
the company must make a capacity expansion decision. 
That is, it must decide at what time (or times) and by 
what amount (or amounts) to invest in capital equipment 
so that production can satisfy demand. The intuitive 
difficulty of this problem is the fact that a change 
in the capacity in one period affects the production- 
inventory schedule for all future periods, thus lending 
combinatorial aspects to the capacity expansion problem. 
Further difficulty arises when capacity must be added 
in indivisible units, each capable of producing large 
numbers of units of output as is usually the case in 
industry. 

Capital equipment itself has many aspects. 


Deterioration might be a significant factor with the 





lThe existence of capacity in excess of demand 
does not preclude a decision to expand capacity since 
there could exist circumstances where such a move 


would lead to lower costs. 





type of capital installed, in which case the capacity 
expansion problem is further complicated by the need 

for equipment replacement. On the other hand, fixed 
capital might have a relatively long life, making per- 
iodic maintenance the only significant factor xc 

the time-span of interest. It is also possible that the 
fixed capital equipment can be salvaged or scrapped if 
warranted by changes in the demand situation. The first 
case can only be approximated by the models which 

follow and then only by allowing maintenance costs to 
reflect the increased investment over time needed to 
maintain equipment at its original capacity. In add- 
ition to not obtaining equipment replacement decisions 
from the models, the approximation is faulty in that 

the increased maintenance costs are applied uniformly 

to equipment in each period, regardless of age. The 
latter two cases are modeled exactly, however, the sec- 
ond by the explicit structure of the maintenance costs 
and the last implicitly by the relevant mathematical 


programming model. 


I.3 Economic Background 


With respect to the economic geography, the 
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problem lies on the boundary of capital theory and the 
classical theory of production. A good discussion of 
the problems arising from the merger of these two dis- 
ciplines is contained in the unpublished Ph. D.thesis 
of Carroll (Carroll 68). 

Apropos of capital theory, a decision must 
be made as to the proper points in time, and the proper 
amounts, for investment in capacity; however, there is 
no requirement present for alternative investments. 
In general, if expansion of capacity in a number of 
different areas were being considered in each period 
with a limited budget, then the extensive work in cap- 
ital budgeting would be relevent.? Further capital 
budgeting problems typically abstract from the micro 
problem of optimizing the form of each investment pro- 
ject, being concerned mainly with the macro allocation 
of funds to projects. The micro-optimization is an im- 
portant part of the work that follows. Extensions of 
the present problem to the capital budgeting sphere 


will be discussed in the final chapter. 


E—————————— MAMMA AIC C - c 


en classical application of mathematical program- 
ming to the capital budgeting problem is given by 


Weingartner (Weingartner 63). 





T. 

In the theory of production as discussed by 
Carroll, consideration of production, inventory, and 
capacity expansion places the problem in between the 
Short term, where capacity is fixed and decisions are 
only made on levels of production, and the long-run de- 
cision of optimal overall capacity level, which ab- 
stracts from all short-run decisions. In effect, the 
mechanism of inventorying production functions as the 
mediator between the short-run minimum cost production 
problems and the long-run minimum cost capacity expan- 
sion problems, serving to "smooth" the extremes in the 
optimal solutions to the individual problems. 

A more satisfying concept of the spectrum of 
decision time-spans is given by Maxwell (Maxwell 65). 
He considers the very short-run as the time-span in 
which all factors of production are fixed in both ac- 
quisition and use. Successive decision time-spans are 
encountered as factors become variable in use and ac- 
quisition, giving a succession of short runs of in- 
creasing length. The long-run is finally achieved when 
all factors become variable in both use and acquisition. 
Against this spectrum, the problem at hand ranges from 
a succession of relatively long short-runs when only 
capital stock is fixed in acquisition and all factors 


are variable in use, to a long-run over which the 





D. 
capital stock 1s variable in acquisition at a succession 
of discrete points. Correspondence with classical 
theory is not exact though, since capital is durable 
and not perfectly flexible as required in the classical 
long-run. The necessity for providing inter-period rela- 
tlonships through the mechanism of inventorying in the 


case of inflexible capital was pointed out by Carroll: 


eeeDO long as the theory of production ab- 
stracted from the particular problems of fixed 
in-puts, It was essentially a mono-period 
theory in which inter-period relationships 
could largely be neglected. Even the long- 
run analyses tended to treat the question of 
determining the optimal capital stock as a 
once-and-for-all decision. Since capital was 
generally assumed to be completely flexible 
in the long-run, questions of the impact of 
one period's capital accumulation decision 
upon that of another were not raised, But 
once capital came to be viewed as a durable 
good, which was not perfectly flexible, 
inter-period (1.e., between long-run periods) 
relationships were a matter of interest. 
(1968, p. 5) 


In the following work, a salvage value of capital 
goods is introduced and constrained to be less than the 
acquisition cost, thus ensuring at least a mild flex- 
ibility in capital. 

Maxwell also presents a variation on the con- 
ventional theory of the firm which is relevant here. 


The contrast is in the derivation of the production 





Un 
function which, in the classical theory, is obtained 


from the combination of varying inputs of consumptive 
factors of production with the fixed factor to obtain 

a graph of maximum outputs over the range of the con- 
sumptive factors. The underlying assumption is, of 
course, that the fixed factor is indivisible in use and 
must be used to the available limit. If instead, the 
fixed factor is viewed as indivisible in acquisition, 
but divisible in use, a more useful approach to the 
firm's production function can be derived. 

Assume constant returns to scale in the pro- 
duction function. Then, for each level of input of 
consumptive factors, the fixed factor is used up to its 
zero level of marginal productivity giving a linear ex- 
pansion path for the firm up to the capacity of the 
fixed factor. Thus, under constant returns, the scale 
of use of the fixed and consumptive factors of produc- 
tion is varied linearly up to the capacity of the fix- 
ed factor, while the proportions of both remain cons- 
tant. Beyond the capacity of the fixed factor the pro- 
portion of the consumptive factor is increased, giving 
returns to proportion in production as opposed to the 
returns to scale obtained in the former range. The 
importance of this concept to the theory of the firm is 


pointed out by Maxwell: 
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Many empirical cost studies indicate that the 
marginal costs tend to be constant very nearly 
up to capacity, suggesting that this range in 
which less than all of the stock of fixed- 
factor services 1s being utilized may indeed 
be a significant portion of the total range 

of outputs that the firm is capable of pro- 
ducing in the short-run. (1965, p. 106) 


For the problem of interest, the range of the 
production function reflecting returns to proportion is 
of no interest since one of the decision variables for 
each period is the optimal level of fixed capital for 
that period. Thus, concern can be focused solely on the 
constant returns to scale portion of the production 
function assumed present up to the capacity of the fix- 
ed factor. 

In the problem which follows, the bundle of 
consumptive factors of production 94 which is needed 
in period i for one unit of output is assumed known. 
Further, it is assumed that the fixed factor of pro- 
duction is used at its zero marginal productivity level 
for this bundle 9. Or any multiple of 1t. The capacity 
W, of the fixed factor in any period can then also be 
expressed in terms of the maximum number of bundles Mn 
which can be used efficiently with the fixed factor in 
the period (in this case the period includes only reg- 


ular time work by one shift). The level of production 





Td 


in period i which is x, must thus satisfy 
x, = W, i = 1,...,N 


These assumptions, when taken in the context 
of fixed capital investment, suppress the problem of 
labor force. That is, although the labor force is a 
non-consumptive factor of production, the necessary 
assumption here is that labor force is always adjusted 
to the exact size required to run the capital equipment 
at the predetermined optimal level x; in period i. An 
interpretation of the problem of capacity in terms of 
labor force rather than capital equipment will be 
given subsequently. 

The objective is to minimize the present 
value of the costs incurred over N periods to meet a 
predetermined demand schedule. It should be noticed 
that the problem is abstracted from pricing consider- 
ations and from determination of optimum output 
schedules; in fact, concern is with the optimum method 
of meeting any demand schedule. To achieve the pre- 
sent value criterion, it is only necessary that the 
individual costs in the objective function of the model 
be discounted by the appropriate cost of capital Py in 


period i. Thus discounting does not have to be modeled 





Were 
but can be introduced explicitly in the relevant model 
coefficients. The subject of the proper choice of the 
discount factor p is not germane to this thesis and 
will not be discussed. 

The major stumbling block to the merger of 
the capital and production theories has historically 
been the problem of fixed capital as a non-consumptive 
factor of production. This question is also lucidly 
discussed by Carroll. To summarize, some method must 
be found to measure the value of a factor of production, 
i.e, fixed capital, which is not consumed in the pro- 
ductive process, but which must be present for produc- 
tion to take place. In this thesis, the value measure 
will be the reduction in cost obtained in the produc- 
tion of a given demand schedule, a measure which is im 
plicit in mathematical programming optimization. 

Finally in the economic context, mention 
should be made of the historic bugbear of marginal 
analysis, indivisibility in acquisition of capital 
stock. This discreteness in acquisition possibilities 
destroys the differentiability of most useful functions 
at best, and their continuity itself at worst. As will 
be described later, a simple adjustment in the integer 


programming formulation will handle this problem 
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fairly easily and, in some cases, even improve the al- 


gorithms. 
I.4 Problem Formulation 


The interval of time under consideration is 
divided into n periods of equal lengths with the first 
period designated as 1 and the nth, or horizon, period 
designated N. It is assumed that demand in period i, 
d,;, is known deterministically for 1 = 1,...,N. Then 
using the following variables: 

X4- production in period i, 1 = 1,...,N 

v47 inventory from period i to period i + 1 

W- capacity added in period i, i = 1,...,N 
the following coefficients: 

p,- cost of a unit bundle 9,0f consumptive 

factors of production in period i, 

DX. d 

h4- unit inventory cost from period i to 1 + 1 
c4- construction cost of a unit of capacity 

in period i 

m, - maintenance cost of a unit of capacity in 

period i 

where coefficients are assumed integer and 


> 
Py Z 0, h, Ž 0, c} aonn aN 





14, 
and the following explicit integer variables: 


Wo- initial capacity (wo 2 0) 


v = initial inventory (vo 250) 


O 
if all functions are assumed linear, the capacity ex- 
pansion-production-inventory problem (CPI) can be 


written as the following integer programming problem: 


N i 
Pi: min Z, = PP ivi teens (E, witwo ] 


ESL. 
TU d 
(P1. 1) X1 - jon Wo ie a... 5N 
XM ec v s di VO 
(P1.2) X4 + V4 L - Vi = d4 i = Dinge le 


where the right-hand side of the constraints is assumed 
integer. Constraints (P1.1) restrict production to be 
less than or equal to capacity in any period, and con- 
straints (Pl.2) require demand to be met exactly out 

of the production and inventory in each period. All 
costs are non-stationary and are explicitly discounted 
to the first period by the appropriate discount rate py, 


which is external to the model as discussed in section 





15% 

Te. 

Initial inventory, v,, and initial capacity, 
Wo, are set explicitly to the values relevant to the 
problem being solved. Initial inventory and initial 
capacity are rather straightforward in the problem con- 
text and no difficulty in obtaining them should occur 
in the real world. Final inventory, however, is more 
of a subjective variable. If vy is allowed to be an 
implicit variable, then Vy ” 0, since if Vy ? O then 
the level of production in periods N, N-1, ... could 
be reduced until VN OU 0, giving a lower cost while 
still minimizing cost (recall that cost coefficients 
are non-negative). On the other hand, Vy can be use- 
ful as an explícit variable as shown by an example. 
Assume, for N = 3, that the demand schedule ís 
df = (1,7,20).3 In this case, vy could subjectively 
be set to some value greater than zero as a surrogate 
for an expectation of continued increase in demand be- 
yond the horizon. 


The objective function in Pl can be re- 


A E IES LCL ALGAE EIEIO LOT SO que AA EE 


3In the remainder of this thesis, an untransposed 


vector is assumed to be a column vector. 
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written as 
min Z4 = i [p4 x4*th4 v4*(c,* Y m,)w, ] * wo: 1 m 
j=] jai J j=) i 


It can be clearly seen from this aggregated function, 
that the total cost of an incremental unit of capacity 
wi is equal to its acquisition cost in period i plus 
the discounted sum of the maintenance costs for one 
unit of capacity up to the horizon. Maintenance costs 
themselves are period dependent ( i.e., independent of 
the age of the installed capacity). 

The CPI model can be extended easily toa 
capacity increase-decrease, production-inventory model 
(CIDPI) as shown below. Let: 

r;~ capacity scrapped in period i 
fi- salvage value per unit of capacity 

scrapped in period i 
all other variables and coefficients are as before and 
assume 


f < C1 j = LN 


i 


N 
L m, < Cy + 


ict 
J kaj 


eq y J > 1 1 — y 
Then the CIDPI can be written as the following integer 


programming problem: 





T7. 


N 


S.t. 
i < 
Na Y 
(P2.2) X1 T Viel = Vi Et di i = ee, Net 
XN * VN.] = dy + VN 
1 > 
X, V, W, F20 and INTEGER 


where the right-hand sides of the constraints are in- 
teger as before, and constraints (P2.1) and (P2.2) have 
the same interpretation as previously. Constraints 
(P2.3) ensure that current capacity is always non- 
negative. 

It should be noted here that P2 contains 
the flexibility of allowing fy < 0 as might be the case 
if, instead of a salvage value for equipment, there 
were a positive cost of scrapping capacity in some 
period (in which case scrapping might be a means of 
avoiding unnecessary maintenance cost on unneeded 
equipment). For fi = 0, P2 is bounded because of the 


initial assumptions made on c, and fs, and for fy < O, 
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P2 is bounded because of the form of the objective 
function. 
As with Pl , the objective function of P2 


can be aggregated as follows: 


min Zo = PROCEDE È mw - Ge E mprg] 
+w’ : In4 
9 1-1 
Pl becomes a special case of P2 if the 

scrapping cost is set to infinity (i.e., f4= -@vyi) in 

P2 . Then r - O0 always and constraints (P2.3) become 
redundant, giving back the form of P1. Two distinct 
models will be retained, however, for both economic and 
theoretical reasons. 

Finally, ín the CTDET model, tne value of 
fixed capital is central to the minimum cost solution. 
During periods 1 to N, this value is reflected in the 
salvage coefficient f. To be complete, though, the 
model must include the value of fixed capital after 
period N, and here two cases are relevant; either the 
firm expects to go out of business after period N or it 
expects to continue production, but its planning horizon 
is limited by dependable demand projections. In both 


cases, a dummy period N + 1 can be added to the model 





19. 
with dy4, = 0. Then, in the former case, fN+1 is just 
the relevant salvage value of the fixed capital, while 
in the latter case, ÍN43 can be set to equal some dis- 
counted estimate of the future value of production 
possible with the installed capacity. 

As for the CPI model, nothing need be done 
since the salvage value of capital has no bearing on 
the problem. That is, capital only has value in the 
problem context to the extent that it can meet the de- 


mands. 


I.5 Related Literature 


A large literature exists in the area of pro- 
duction planning and there is at least a modest one on 
capacity expansion. But the combined problem, though 
economically significant, has been given relatively 
little attention, and no attention at all has been 
given to the necessity of providing good integer solu- 
tions to such a problem. A related literature closely 
allied with the combined problem does exist, however, 
in the area of employment, or work-force, smoothing. 
Each of these areas will be discussed subsequently. 

Mention must first be made of the relation- 


ship among the existing capacity expansion literature 
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economics, and reality. Any demand-related problem 
naturally depends on a reasonable projection of demand, 
either stochastic or deterministic, within a believable 
horizon. Since forecasting is, at best, a risky busi- 
ness, the veracity of any demand projection must have 
an inverse relationship to the length of the horizon 
over which the projection is made. A reasonable figure 
for most companies would probably not exceed five years. 
Within such a time span, however, costs of production 
and inventory become economically significant, while 
beyond this span, any quantitative approach rapidly 
approaches speculation. The only valid conclusion then 
would seem to be that a realistic, quantitative 
approach to any capacity problem in the production sec- 
tor must include the effects of the production and 
inventory policies. 

A further economic and realistic condider- 
ation in the capacity expansion problem must be the 
fact that capacity, in whatever form it may take, is 
usually obtainable only in fixed increments. Further- 
more, an incremental acquisition of capacity could, by 
the nature of the machine or whatever that is added, 
correspond to a large incremental capability in terms 
of units of output. Thus, although the roundoff of an 


exact, non-integer solution may give an answer that 





21 
is "clese" to optimal in terms of the demand profile, 
the error in terms of cost could be very large. It 
must be concluded then that, at least in the context of 
capacity expansion, integer solutions are important. 

Perhaps the classic early delineation of the 
capacity expansion problem was presented by Manne 
(Manne 61) in his analysis of the basic problem for 
both deterministic and stochastic demands. Early 
linear programming models of the problem were pre- 
sented by Fetter (Fetter 61) and Veinott and Wagner 
(Veinott-Wagner 62). An equivalent form of the 
latters' model and a corresponding network formulation 
can easily be obtained from P2 by setting production 
and inventory costs to zero. Kalro and Arora (Kalro- 
Arora 70) extended the basic linear programming form- 
ulation to handle stochastic demands. 

A dynamic programming approach to capacity 
expansion is given by Howard and Nemhauser (Howard- 
Nemhauser 68) for deterministic demands and extended 
to the stochastic demand case by Montgomery (Montgomery 
71). An extension of the dynamic programming model 
for deterministic demands to handle concave costs of 
expansion is made by Manne and Veinott (Manne- Veinott 
67). None of the above work takes into account the 


production-inventory problem, however, and thus will 
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not be considered as directly antecedent to this thesis. 

Zabel (Zabel 63) was one of the early authors 
to look at the total problem, but his assumption of no 
inventory effectively reduces the problem he solves to 
one of capacity expansion. Hinomoto (Hinomoto 65) 
essentially retained Zabel's no-inventory assumption 
and extended his work to include technological impro- 
vement, thouch he couches the work in terms of replace- 
ment theory (failing to mention Zabel). 

Iglehart (Iglehart 65) seems to have analyzed 
the first complete model for the CPI problem. He ex- 
tends Zabel's work to allow inventory and restricted 
back-ordering using analytical techniques to establish 
optimal capital accumulation and production policies. 
He poses what seems to be a serious restriction, how- 
ever, when he assumes that the consumptive factors of 
production are always spread evenly over all fixed cap- 
tal, Mn effect eliminating the possibility of idle 
capacity. Carroll (Carroll 68) extended Iglehart's 
work to stochastic demands and pricing considerations 
using dynamic and non-linear programming techniques to 
establish optimal capital investment and production 
policies, 

The most significant general formulation of 


the CPI problem to date, from a quantitative viewpoint, 
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has been the work of Cobian and Lele (Cobian-Lele 70). 
They give a general dynamic programming formulation of 
a multi-plant, multi-inventory facility firm faced 
with deterministic demands and show how the model can 
be hierarchically decomposed to simplify the solution 
process. Unfortunately, even when only one plant and 
one inventory facility are assumed, the model still 
contains two state variables and gives non-integer solu- 
tions. 

The area of recent work which appears most 
relevant to the area of capacity expansion seems to 
have been done in the related literature on employment 
smoothing.  Interpretationally, work force level can 
be viewed as a measure of capacity; the assumption 
made in this literature being that the fixed capacity 
level is given in each period, as opposed to assuming 
that the work force itself is optimally adJusted to 
the level of output, as mentioned in the previous sec- 
tlon on economics. The greatest objection overcome in 
the employment smoothing literature is that, by nature 
of the problem, the production-inventory policy is 
always included. To date, however, integrality has not 
been considered an important aspect of the problem, 

An excellent delineation of the problem and methods of 


solution plus a bibliography through 1967 is given by 
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Silvenw4Silver 67). 

A linear programming model which is close in 
form to P2 and a heuristic for "rounding-off" the 
linear programming solution are given by O'Malley, 
Elmaghraby, and Jeske (O'Malley et al. 66). As is well 
known, and demonstrated by example later in this thesis, 
rounding-off can be far from optimal and a very costly 
approximation to the optimal integer solution. 

The major difference between P2 and the 
model of O'Malley et al. is the non-stationarity of all 
costs in the former, as opposed to the complete station- 
arity of costs in the latter. In addition, the capa- 
city decrease coefficient in P2 can be positive or 
negative, within the bounds of the preconditions, to 
correspond to salvage value or scrapping cost re- 
spectively, while the corresponding term in the O'Malley 
model is a firing cost and constrained to be positive. 
Lastly, P2 is an integer program as opposed to the 
linear program of O'Malley. Still, O'Malley's model is 
not a special case of P2 since it has a term for over- 
time production. This, however, is a simple extension 
of P2 done in Chapter III, at which time the O'Malley 
model will become a special case (except for integral- 
ity). 


The balance of the recent major work in 
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employment smoothing has been done in two companion 
papers by Lippman, Rolfe, Wagner, and Yuan (Lippman et 
al. 67a and 67b) and a later paper by Lippman and 
Yuan (Lippman-Yuan 69). The model addressed in the 1967 
papers is a generalization of the linear model of 
O'Malley et al. with convex-like, non-stationary pro- 
duction costs, increasing and non-stationary inventory 
costs, and v-shaped, stationary smoothing costs. De- 
mand-dependent upper and lower bounds on work force 
size are given &nd optimal policies are developed for 
the special cases of monotone increasing and monotone 
decreasing demands. Algorithms are then given for the 
special cases of monotonicity and the added restriction 
that all costs are linear, though non-stationarity is 
retained for production and inventory costs. 

In the 1969 paper, similar analytical work 
is done on a model which is the same as that of 1967 
with the exception that smoothing cost is permitted to 
be non-stationary. This paper also points out that the 
upper bounds and algorithms developed in the 1967 
papers are invalidated under the more generalized 
smoothing costs. No further algorithms are presented. 

The models in this thesis are special cases 
of the work of Lippman et al. only insofar as the more 


general form of the latter's production and inventory 
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cost functions are concerned. Integrality and the 
existence of algorithms in this work seem to make 
further comparison unnecessary. The most cogent point 
in the above discussion would seem to be the recogni- 
tion of the interrelationship between capacity expan- 
sion literature and employment smoothing literature, 


about which more will be said in the final chapter. 


I.6 Plan of the Work 


In Chapter II, Pl and P2 will be trans- 
formed to the equivalent problems on which the algor- 
ithms given in the chapter are based. Bounds are de- 
veloped here for the variable which acts as a surrogate 
for capacity. Finally, more efficient bounds and 
algorithms are presented for special monotone demand 
and cost cases. 

Chapter III presents a network formulation 
of the entire problem for the special case when cap- 
acity increase and decrease costs are equal in each 
period, but are allowed to be non-stationary. The 
basic models are then extended to include indivisibil- 
ity in acquisition of capacity and piecewise-linear, 
convex production costs that provide for consideration 


of multi-shift work and overtime. 
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A discussion of similar integer programming 
algorithms and a reinterpretation of the results of 
this thesis for the employment smoothing problem are 
presented in Chapter IV. Also presented here are the 
results of test problems run with a FORTRAN IV coding 


of the algorithm along with a summary and conclusions. 





Chapter II 


I| -Introduetion 


Models Pl and P?2 of Chapter I could be re- 
formulated to eliminate the inventory variable v, in 
which case Pl would have e2N constraints and 2N variables 
and P2 would have 3N constraints and 3N variables. In 
the capacity expansion context with a period of one 
year and N = 5, neither problem would be very large and 
existing all integer algorithms could be used. On the 
other hand, in the employment smoothing context with a 
period of one month and N = 24, where furthermore the 
CIDPI model is most relevant, the problem size is out 
of the range of all except branch and bound type al- 
gorithms. 

Use of implicit enumeration seems out of the 
question since the variables of Pl and P2 can take on 
rather large values, and expansion to O0 - l variables 
would quickly exceed the capacity of existing algorithms 
for even small values of N. Primal and dual cutting- 
plane algorithms have the disadvantages of highly un- 
predictable convergence and, for the latter, lack of a 
current feasible solution until optimality is reached. 


Although there is the possibility that one of these 
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other methods possesses an extremely efficient adapta- 
tion for this special problem, this thesis investigated 
the most obvious approach- branch and bound. In this 
chapter, special branch and bound algorithms will be 
derc oped tor Pl and P2 which exploit the structure in- 


herent in these problems. 


In section II.2, a linear transformation is 
done on Pl to obtain an equivalent form P3 for which an 
algorithm can be more easily developed. The nature of 
the algorithm is discussed subsequent to the transform- 
ation. When treated as a function of certain of its 
variables, the optimal value of the objective function 
of P3 has convexity properties over certain ranges of 
those variables which are established in section II.3 
by a series of lemmas, theorems and corollaries. The 
properties of the integer program P3 plus the results 
from section II.3 lead to the set of bounds and op- 
timality tests given in section II. which are sub- 
sequently used in the branch and bound algorithm of 
section II.5. 

A similar development is followed for the 
CIDPI model P2 in sections II.7 and II.8. The linear 
transformation used is the same and the resulting 
model P5 has properties similar to those of P3 both in 


problem structure and objective function convexity. 
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The modifications to the basic algorithm of section 
II.5 needed for P5 are given at the end of section II.8. 
For both transformed models there are special 
cases of demand and cost structures which lead to major 
simplifications in the algorithms. These structures, 
their properties, and the attendent modifications of 
the algorithms are given in section II.6 for the CPI 


model and section I1.9 for the CIDPI model. 


II.2 Transformation of Pl 

A set of necessary and sufficient conditions 
for the feasibility of a W can be obtained from con- 
straint sets (Pl.1) and (P1.2). The method in general 
is to solve the inequalities (P1.1) for x by adding 
the appropriate terms to each side. Then for any 
€? ....N, sum the first k constraints of (P1.2) to 


obtain a constraint of the form 


and substitute for the x, from (P1.1) giving an in- 


equality of the form 


k 
Eid; - Vo - kw 


> 





Bale. 


Repeating for each value of k gives the set 


z 
= 
i 
< 
e 
< 
t 
z 


ho 
X 
ps 
+ 
€ 
ho 
i 
< 
ho 
i 
© 
+ 
© 


Nw, + eee t Wy a VN - d) ł eee + dy hr Vo - Nw, 


If a vector of surplus variables t 2 0 is 
added to (2.1) and the v, and t, in each constraint of 
(2.1) are aggregated as 

PE V4 TU 


di i = Je M 


then y 2 0 since V 2 0 and t 2 0, and (2.1) can be 


written in vector form as 


where 
l 0 0 2c On 
c L 0 DES ie 
j à 1 o e e 0 0 
B = 0 o $ e 9 e e e 


N (N-1) (N-2) . . . 2 1 


and the form of the vectors is obvious. 
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Solving for w from (2.2) gives 


^ 


NM UI A es ie 
w * B^^d * B^^y - B-lV, - B-lw, 


where 
teen O20 4 «= 0 0 0 
2120 0 2-7. 9240 O 
um? wm. ..0 O0 O0 
0 1-21...0 00 
BI = MENU uU... 


Expansion gives the following set of relations 


w, = d, - dy t y2 -= ey, + V. 


2", 3) 
w3 = d; = dp + y} - 2y5 * y4 
AS YN t XN-5 


which expresses w explicitly in terms of demand, 
initial inventory, initial capacity, and the newly in- 
troduced aggregate variables y- The structure of the 
right-hand side is rather unique and, as will be seen 
shortly, extremely useful. The above transformation 


will be referred to as Standard Transformation l. 





33% 

An interesting economic interpretation can be 

given to the variable y. If overcapacity is defined as 
the cumulative excess of capacity over cumulative de- 

mand up to and including period i, then Ys is the over- 

capacity in period i. For instance, if N = 3, a! = (3, 


5,6), and wi = (5,0,0), and MON VASO then y, & y». v 2 


O 


and o l. Further, in any period 1, v, is the por- 


i 
tion of overcapacity Y 4 which is used to produce for 
demand after period i and stocked as inventory, and ty 
is the cumulative amount of idle overcapacity as of the 
end of period 1. Thus, if VN ” 0, then yy is the total 
idle capacity in the model. 

Model Pl can be transformed by substituting 


ror Ws from Standard Transformation 1 to give the 


following equivalent model.! 


P3: min Z4 = N max + faves + ety, | + K 
set, 
< d. 4 
a A ere YG 
(P3.1) 2 
T dy t Y4 ~- Yil diet Be oe lc ey N 





lPor the remainder of this thesis, the objective 


function of P4 will be referred to as Zj. 
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(P3.2) X4 + rem m MET - di i = estesa 


iv 


Ji Wo + Vo => di 
> 
2 i 2 d d = 
Vs T ye) Yiz2 cmm dedi = 4 i SN 
X, V, y 20 and INTEGER 


where constraint set (P3.3) comes from substitution in 

w — 0 in P1; the restriction of y to integers comes from 
the restriction of w to integers in Pl and the fact 

that the right-hand side of constraint set (P3.3) is in- 


teger by assumption, and where 


' = = = = = 
E Cy 2€141 + Cayo + m, Ms 44 i de 


ES =D Sen VEN) T 


rm 


Nel 
K = E, LC 047014 1*4) dy I+ Coy tmy) dy-eywotlca-e1-m1 IV, 


If, in P3, a feasible y is given then the problem 





Dar 


becomes 
y N 
SU 
« 
= + = 
x. d, yy NO 
(P4, 1) e 
Ly x d4 t Yi — XE i = zoe | 
Xi - Ya = qd, - Va 
(P4,2) Xi + yisi = E = d, i = 2p... ,N-l 


Xy * Vy-1 7 Gy t Yn 


2 0 and INTEGER 


x, V 
which Just is the production-inventory problem, para- 
metric in y; but which, more importantly, can be shown 
to be a transportation problem as was first pointed out 
by Bowman (Bowman 56). Specifically, let 

X147 production in period 1 for demand 
> 


in period j, jJ- i 


then P4 can be written as 


N N 
PY! + min Z' = E Eu a,,X 
ho i91 jei 13 1J 
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Se G's 
qd, +t y, = v Je v1 
(Ph 1") LI A uu ou 
jzi Ay AT j = ETT PN 
J dj - Vo Jeu. 
o. 2') 5 Xij = d a = se]. 
i=l AN + VN J = N 
x 20 and INTEGER 
where 
j-1 
as = Pi + En J n een 
«11. — Pg 


and x is integer by the properties of the transporta- 
tion problem. 

The impact of this development is that, in a 
branch and bound algorithm enumeration would have to 
be done on only a subset of the problem variables, i.e. 
y. For each feasible y generated, corresponding to a 
node in the enumeration, a transportation problem can 
be solved. This, in effect, implicitly fathoms each time 


what would have been a significant portion of the 





©The lower summation range in (P4.1') could be ex- 
tended to 1 and the upper range in (P4.2') to N if 


814 =00,Vj < i, i = 2,...,N. 
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enumeration tree had branch and bound been tried on 
the entire problem. Furthermore, the number of trans- 
portation problems solved can be cut to a small fraction 
of the number of nodes actually visited by using suit- 
able tests and bounds. In effect, the branch and bound 
is a search for the minimum cost combination of a 
feasible y and its corresponding optimal production- 
inventory schedule. 

Lower bounds for y can be obtained directly 


from constraint set (P3.3) which can be solved for y as 


iv 


iv 


1 2 O 


This set of lower bounds on y is integer since all 
terms on the right-hand side are integer. 

Upper bounds on the problem variables must 
exist since the model is bounded . The following lemma 
will give a set of global upper bounds on Y. 


Lemma 2.1 Global upper bounds for y are 


N J 
(2.5) y, ^ J- E d, - Ey IVY DI dan 
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Proof: Eor) œr from (2,2) 


O O 


and by inspection the right-hand side is maximized when 
Wa takes on its maximum value. The maximum capacity in 
period 1 would be a level sufficient to produce all 


future demands in period 1. This value is 


= 
p 
t2 
B 
{ 
< 
i 
= 


1-1 P 9 
and so 
max y, = NUT - dà 
or 
EN 
cac 


In general let J = k, ke{2,3,...,N} 
O O 


k 
Yy = kw] + (k-1)w, +... + Wy - 2,41 + + kw 


Again the right-hand side is maximized when Wj 18, and 


from the argument above 


N k 
EB INS a (k=1)v 
K 12] 1 j=] 1 O 


so the result is true in general. 
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As wlll be discussed later in this chapter, 


upper bounds are not really crucial to the algorithm 
since more powerful termination rules will be employed. 
A special case which has an optimal solution with 


m O0 will also be investigated subsequently. 


IAS Properties of P3 


Since a branch and bound algorithm is to be 
developed it is important to have efficient bounds and 
termination rules in order to perform as little enumer- 
ation as possible. The following lemmas will establish 
useful properties of some variables of P3. 

Lemma 2.2 If all y, are fixed for i = 1,...,N-1 and YN ds 


lneremented by n to 


YN mE 9 EM o 
then the new value of Wy is 
l = 
Wy Wy tn 


Proof: From (2.3) 
Tw = 
Dx Ee - duy 


(yy) ~ 2¥y 4 * Yu o * Ay > A 


w rn 


N QED 


For Lemmas 2.3 and 2.4, Corollary 2.2.2 and Theorems 


2.1 and 2.2, when given a set of Yi» ie eer duis 
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assumed that the remaining Y 3k are at their lower bounds for 


Lemma 2.3 If for given Y4 5 rss ees, the lower 
bounds of Y 3+k are positive, k = 1,...,N-J then Wiek * O, 
kradl s.. N=]. 

Proof: For any keí1,2,...,N-J) from (2.4) and the 


hypothesis 


> 
Yjek 7 ?Yj«k-1 - Yj«k-2 * dj+k-1 7 đj+k > 0 


Further, if Yj+k is at its lower bound, equality holds 


on the left. The expression for Wirk from (2.3) is 


"uL "LYj3«k 7 é3j)4k-1 * Yj4k-2 +t dj+k-1 7 dj+k 


but the right-hand side has just been shown to be zero. 
Thus 
Wi+k * 0 
QED 
Lemma 2.4 If for given Y,» 1 = 1,...,J, the lower 
bounds of yj4y are positive, k - 1,...,N-J, and y: = 


y 3 * n; nei... then 


w! xz Wy +n 
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Sud Tor 1200. ¿N=] , the lower bound on Y ek is 


Y 3+k + (k+1)n 


and 


Proof: From (2.3) 


a = yj - ?Y 421 + y + an - di-1 
SO 


a = (y , +n) - 2Y 3-1] + Yj-2 + d, - diy a Wy +n 


which proves the first assertion. The second and third 


are proved by induction on k. For k = 1, from (2.4) 


> 
n - 2(y,*n) - yyy + dj - 4541 = Yj+1 + 2n 


and from (2.3) and Lemma 2.3 


' = *2n) - 2(y,*n) * +d -d 
"Pn rng a» pru eg J+1 j 


Assume the relationship holds for a general ke{2,3,..:.; 


N-j-1}. Then for k+l from (2.4) 
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* 2ly j4* (k*1)n] - LY 44 k-1*tkn] * ayy, - d 


ut 
3+k+1 3+k+1 


N 


Yj+k+] + (kt2)n 


The final assertion is then true by Lemma 2.3. 
QED 
Remark: From the relationships (2.3) and (2.4) it can 
be seen that a change in Y 4 leaves Wises a 4] and Y1» 
0,9 j 1 unaffected. 
Now let 
Cly y) = Capacity cost as a function of y; 
PI(yj) = production-inventory cost as a 
function of yj 

Ty 5) = total cost as a function of y, 

so 


T(y 4) = Cly;) + PI(y,) 


Then the preceeding lemmas can be used to establish 

the following cost properties for P3: 

Theorem 2.l If for given Ys. il. seg eae lower 
bounds for Ysa, are positive k = 1,...,N-j then C(y ,*n) 
is linear and non-decreasing in n for ne{1,2,3...}. 


Proof: From Lemma 2.4 for y. = y,*n, rhe fl 253.065 


i: = Ws +n 


Wir = 0 Fame eon Nad 
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mmus from Pl 


A A ACTES S E O APN: 
E - C y y m y) eee Si ee = 


J 


QED 
Before stating the next theorem, a well 

known result for the linear program associated with 
P4' must be mentioned. Since the linear program 
satisfies the sufficient conditions for unimodularity, 
then all extreme points of the convex hull are integer 
(Hu 69). This allows use of the properties of the 
general linear program in the following proof. 
Merci nie). + and for any Jell,2,...,N) 
Lf y. - Y 4 + n and for each value of n the y, are set 
at their lower bound, i = jtl,...,N, then ey is 
convex inn and the Yi A A A 

Proof: Assume the following linear programming 


problem 


and the two parameterizations of b 
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o 


(p2) Ax ^ b + (1-A)E 


ol 
+ 
mi 


(p3) AX = 


< 
Were 0 = \ = l,€ is any vector which maintains feasi- 


bility in (p3) and the optimal solution of (pl) is Xi 


and of (p3) 1s X3 where 


Te = 
C Xi f (b) 


P E f (b4£) 


are the optimal objective function values as functions 


of the right-hand sides. Now let 


i 
PA 


tS tom AX. + (1-A)x 
( ) x, X (1 ) x3 0 
and since the problems are linear it must be that 


(2-77) ey, = cl, + (1-A) e7X, 


Multiplying the constraints of (pl) and (p3) by the 


appropriate constants gives 


AX, = Ab 


ÍA 
ÍA 


(1-A)AxX_ = (1-XA)(b+E) 


3 
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and since à is a scalar, this can be rewritten as 


AAX] = Ab 


A(1-A)Xq = (1-A)(b+E) 
Adding these two constraint sets gives 


A[Ax, + (1-4) x3] 2 b * (1-A)€ 


which, from (2.6), 1s Just 


AX? z b * (1-A)€ 


IA 


IA 


which proves that the convex combination of the op- 


timal solutions to (pl) and (p3) is feasible to the 


convex combination of right-hand sides. 


the optimal solution to (p2) 
~~ < -T- t 
(2.8) Č X, c x) 


since it is a minimization problem. 


From (2.7) and (2.8) 


eIx, $ jcTx, « (i-3)87X 


3 


Now 1f x. 1s 


tA 


2 


tA 
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producing the general convexity result and, by the 
remark preceding this theorem, the desired result for 
PI(y ,+n). 

QED 
The following corollaries follow trivially. 
Corollary 2.2.1 PI(y,tn) is convex for ne{1,2,3,...}. 
Corollary esenec If for given Yas IM. sj. the 
lower bounds of Y 34+k are positive, k = l,... N-J then 


Erin) io convex inn for ne{l,2,3,...}. 


J 
Proof: From Lemma 2.4 the capacities in periods 
J+1 through N are equal to the capacity in period j. 
The result follows from this and Theorem 2.2. QED 
The major result can now be given as follows: 
Theorem 2.3 If for given ZE Dac 0, whe lower 
bounds for Y 4k are positive, k= l,..., ,N-jJ, then 
T(y,*n) is convex in n for ne{i,2,3,...}. 

Proof: The proof follows from Theorem 2.1 and 
Corollary 2.2.2 and the fact that the sum of convex 


functions is convex. QED 


Corollary 2.3.) T(yytn) is convex for ne{1,2,3,...}. 


II.4 Optimality Tests and Bounds for the Enumeration 


The results of section II.3 can now be used 


to develop bounds and optimality tests for the 
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algemithm. The logical first step is to solve the 
linear programming problem for Pl and use 
Optimality Test 1: If the optimal linear programming 
solution is all integer, then this is an optimal solu- 
tonto P3. 


Otherwise let 


A. - optimal objective function value of the 
linear program associated with the integer 
programming problem 

Zip = optimal integer program objective function 


value 


It is a well known fact that 


IA 


Z* Z* 
er Ir 
and this is used as 
Bound 1: A lower bound on total cost is Zip’ 
The algorithm will use Bound 1 in the follow- 
ing optimality test. 


Optimality Test 2: If Zjp is the current integer pro- 


= % 
Z Ip EN 


then Z lp is an optimal solution to the integer 


gramming solution and 
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program. 3 
The next bound uses the minimum production- 

inventory cost, PI*, to obtain an upper bound on 
capacity expansion cost. PI* can be obtained by 
solving P4' with the values of the right-hand sides 
of constraint set (P4.1') equal to infinity. This 
corresponds to putting an infinite supply at the 
sources of the transportation problem. 
Lemma 2.8 If 23 is a current optimal solution in the 


enumeration then 


e > 
N 


Z' - PI* 
3 


is an upper bound on capacity expansion cost. 
Proof: Assume there is a capacity expansion cost 


at such that 


and with a corresponding 2. such that 


Ji 
aS! 
3 3 


But this implies that 


AREA AAA DAA nn nana 


3[a] indicates the next integer greater than a. 
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pr = zt ~ al < 2! ~ a = PI 
3 3 
which is a contradiction of the construction of PI*, 
QED 
Using the above gives: 
Bound 2: If fe is a current optimal solution of the 
integer program, then an upper bound on capacity ex- 


pansion cost a is 
a Í Z! = PI® = & 
3 


Lemma 2.8 can be useful in another way. If 


the current upper bound on capacity cost is a then 


certainly an upper bound on period 1 capacity is“ 


{A 


^ 
a 





+t 
e m 
m l 


since any larger value of Wy would violate Oe A 


similar bound for Wo is 


R N 

cl DT 
Ww s BR ——— 
mee na) 
"ON o i 





MET indicates the greatest integer smaller than a. 
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In general the following bound can be used. 
Bound 3: For the current a, upper bounds on capacity 


expansion in each period are 


I^ 


“7 
cıt 


j» 
MZR? 


UE 


a ir + y ) 
a= c mi )w 
< 1=] 1 x=1 k’ i 
DE m A ANN J N 


N 
SIR Te 

The following bounds use the results of the 
prevlous sections, 
Bound 4: For any feasible values of Yi, 17 1,...,J], 
lower bounds on the values of Yi» l 54-0... Y can be 
obtained from (2.4). 
Bound 5: Global upper bounds on Yas itt (1 tengo can 
be obtained from (2.5). 
Optimality Test 3: If enumeration has been completed 
with NE at its global upper bound or the upper bound 
inferred from Bound 3, whichever is lower, then the 
current optimal solution is optimal for the integer 
program. 

The following bound on yy uses the results 


Olmcorollary 2.3.1. 





DI 


Bound 6: If for-a given set of Woo dtu rl 
< 


then the enumeration on YN for this set of y; can be 
fathomed at (yytl). 

The following bound uses the result of 
Theorem 2.3. 
Bound 7: If for a given Yu» Yo»... >Y3> the lower 


bounds on y, are positive for 1 = jtl,...,N and if 
T(y,) = T(y,*1) 


then the enumeration on y; can be fathomed for this set 
of ZP (= eel at (y¿+1). 

Bounds 8 and 9 make use of Theorem 2.1 and 
Corollary 2.1.1. 
Bound 8: For given Y4» 1 = 1,...,J, if the lower bounds 
Dare pOsative for i= jt1,...,N and if Cly 4) _ a, 
then y; is fathomed for the preceding set of Ys. 
Bound 9: If for given y;, 1 = 1,¥..,N-1, Clyy) 24, 
then no further enumeration on YN is necessary. 

Since the algorithm conducts the search by 
enumeration on y, most often, Bounds 6 and 9 become very 
powerful tools. Bounds 7 and Y are also extremely 


powerful because their conditions will be satisfied 
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often during enumeration on the later period variables 
due to relationships (2.4). The usefulness of these 
bounds diminishes the need for the global upper bounds 
(2.5) since such bounds must include all cases where 


Bounds 6 through 9 would be useful. 


II.5 Algorithm for P3 


step 0 Initialization Set the value of current optimal 
cost, a and the current bound on capacity cost, a, to 
infinity: 

Step 1 Solve the associated linear program. If the 
solution is integer, stop. (Optimality Test 1) Other- 


wise set the linear programming bound to (Bound 1) 


7% 
LP 


Step 2a Solve an uncapacitated transportation problem 


and go to Step 2a. 


to get the minimum production-inventory cost PI* and go 
to Step 2b. 

Step 2b Find global upper bounds for yy, 1 =1,...,N 
from (2.5). -Go+to Step 3. 


Step 3 Start the enumeration by setting 
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y, ^ max E 


and finding lower bounds for the remaining values of yy, 
4 = 2,...,N, from (2.4), using this value of y,. The 
pointer (PNT) is N. Go to Step 5. 

Step 4a PNT = PNT - 1. Go to Step 4b. 

Step 4b With PNT = J, let y, - y,*l and, 1f j « N, find 
new lower bounds for succeeding y, from (2.4). If j Al 
and any lower bound exceeds an upper bound (from Bounds 
3 and 5), yj is fathomed so go to Step 4a with PNT = j. 
If j= l and some lower bound exceeds an upper bound, go 
to Step 10 (Optimality Test 3). Otherwise go to Step 5 
with PNT = j. 

Step 5 With PNT = J determine capacity cost a = cCy D. 
If a > à and either J s N, or J « Nand y, » 0,1 -J41, 
...,N go to Step Ma with PNT = j (Bounds 8 and 9). 
Otherwise if a > A go to Step 4b with PNT = N. If 

geo and si je N with y_ > 0, 10 J+1,...,N go to 

Step 6 with PNT = j (This sets up convexity test for 
Bound 7). Otherwise go to Step 6 with PNT = N (This 
sets up convexity test for Bound 6). 

Step 6 With PNT = j solve a transportation problem 

to obtain the related optimal production-inventory cost 
PI(yj). Then total cost 1s T(yj) = C(y¿) + PI(y,). Go 


to Step 7 with PNT = Jj. 
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BECDOADUWICth PNrae i, if J" Nor Jos N and YI 0, 
134-0] 41,...,N and TCy 4) > T(y,-1) then y is fathomed 
so go to Step 4a with PNT = j (Bounds 6 and 7). If 
Tlyy) = T(y ,71) = D in any case go to Step 9. For all 
other cases if T(yy) < : go to Step 8, otherwise go to 
Step 4b with PNT = N. 
seep 8 If Tly 3) =: ER go to Step 10 (Optimality 
Test 2). Otherwise go to Step 9. 
step 9 Either include T(y4) = Un on the list of 
current optimal solutions or replace the current list 
with 2 = T(y4), whichever is appropriate. Update in- 
termediate bounds from Bound 3 and go to Step 4b with 
PNT = N. 
Step 10 Decode the current optimal solutions and STOP. 

The algorithm, as written, will also find 

all alternative optimal solutions if termination is not 
due to Step 8. The flexibility of choice among solu- 
tions would normally be desirable in a management 
situation. If the alternate optima and en running 
times are not desired then the following changes to the 
algorithm are needed: 
Step 5 With PNT = j determine capacity cost a = Clyjz). 
If a 2 8 and either j = Nor J < N and yy > 0, i = J+1, 
...,N, go to Step 4a with PNT = j. Otherwise ifa? â 


go to Step 4b with PNT =N. Ifa < a and if j < N with 
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y, 7 0, 1 = Jtl,...,N go to Step 6 with PNT = j. 
Otherwise go to Step 6 with PNT = N., 
pep i With eNTe<.), Jf jJ $*N or J < N and y4 ? 0, 
i= j+l,...,N and T(y,) - T(y,-1) then y, 1s fathomed 
so go to Step 4a with PNT * j. If DO < Be go to 
Step 8, otherwise go to Step 4b with PNT = N. 
Step 9 Replace the current optimal solution with 
23 = T(y4). Go to Step 4b with PNT - N. 


II.6 The CPI Model Under Special Demand and Cost 


Assumptions 


In this section two special cases of the CPI 
model will be discussed, one involving monotone non- 
decreasing demands and the other monotone non-increas- 
ing demands. The former leads to significant simpli- 
fications 1n the CPI algorithm, while the latter isa 
trivialwcasewof thes@PI model. 

The following two cost assumptions hold in 
this section: 

(a) The cost assumptions of Pl hold. 
(b) The cost coefficients of Pl (Chapter I) are mono- 
tone non-increasing with time. 


Assumption (b) allows the possibility of discounting. 
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Monotone non-decreasing demand 


The case of demand which is monotone non- 
decreasing is likely to occur often in the context of 
firms interested in solving capacity expansion problems. 
The main result for this case is contained in the 
following theorem. 

Theorem 2.4 If assumptions (a) and (b) hold, demand 
is monotone non-decreasing, and 
< N 
Nw. ven NC 
then there 1s an optimal solution witn yy = 0. 


Proof: Assume in contradiction that 


It can also be assumed, as discussed in Chapter I, that 
ae O. But this just implies that there are a> 0 
idle units of capacity in the model. Consider the first 
unit of idle capacity. One of two cases can be assumed 
to exist, 

Case I- The idle unit was added in period j instead of 


k where it was first used, where j < k SN. 


Case II- The capacity was added in some period previous 
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to ], was used up to and including J-l, and became 
Fale» in j. 
For Case 1 to happen either the cost of adding capacity 


in J must be cheaper than k, i.e. 


N 
n m4 


N 
er + 7 


peo, t 


which, since Cs 2 cy (assumption (b)) implies that 
K-1 


d < 0 contradicting assumption (a); or the cost is 
equal to that in k in which case the idle unit could be 
removed in jJ and added in k with no change in cost. 

For Case II it must be that Wi = 0 since if not, then 
Case I applies. So this unit of capacity must have 
been built in some previous period q and used to pro- 
duce from q through J-l. Now at least one of these 
(j-q) units of production must be inventoried into 
period j since if ey = O then since X1-1 > Xj thas 
implies that dı_ı > d; contradicting the demand assump- 
tion. Furthermore, by the same argument, V4 2 1, 1% 
q,qtl,...,j-l so the idle unit in jJ was used to produce 
for inventory in period q and at least that unit of in- 
ventory was carried into period J. The production 
schedule of periods q*1 to j-1 will not change if the 
unit of capacity is added in q*1 instead of q, however 


the unit which was inventoried from q to j could be 
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produced am jJ instead, eliminating the idle unit of 
capacity. With idle capacity the incremental cost con- 
sidering just the one unit of production in question 


was 
N 
(2.9) cet £ 


while adding the capacity in qtl and producing in j 


gives 
N 
ep + L + 
( ) Sa ee Pj 


Subtracting (2.10) from (2.9) gives 


= > 
(C g-Cq+1) + mg + (Pg-Py ) + E hy = 0 


1 
q 
by assumptions (a) and (b). Thus, in both cases, idle 
capacity can be reduced with, at worst, no reduction in 
cost. 

If idle capacity still exists, the argument 
can be repeated. 

QED 

Corollary 2.4.1: If the assumptions of Theorem 2.4 
hold then all capacity is used in each period. 


The use of Corollary 2.4.1 leads to an imme- 
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diate simplification of the production-inventory prob- 
lem P4. Since the right-hand sides of constraint set 
(P4.1) are Just capacities, equality holds in all N 
constraints and the x, are determined trivially for a 
given Y. Constraint set (P4.2) becomes N equations in 
N unknowns (since v 


is explicit and v,, is either zero 


O N 


implicitly or set explicitly) and is thus deterministic. 
So the subproblem in the enumeration can be solved by 
inspection each time. 

Another observation follows from Theorem 2.4 
and Lemma 2.4.1. If all capacity is always used it can 
certainly never be feasible to build capacity in any 
period sufficient to produce all remaining demand. In 
fact, the upper bounds on capacity increase are much 
smaller than those given by (2.5). Instead of stating 
them explicitly, they can be handled implicitly in the 
algorithm by the following: 


Bound 10: For given y,, 1 =1,...;,j, if the increment 
oo + 
= Y 3 l 


causes the succeeding lower bounds to contain yy ? 0; 
then y 4 is fathomed for the preceding set of y4. 
Bounds 6 through 9 are now useless, since 


their conditions are never met in this case, so 
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Bound 10 takes over as a main fathoming rule in the 
algorithm. Furthermore, the global upper bounds from 
(2.5) are irrelevant to this special case and can be 
ignored. 

The basic algorithm requires the following 
modifications: 
1) Drop Step 2b. 
2) Step 4b* With PNT = jJ let a = y,*l and find new 
lower bounds for the succeeding elements from (2.4). 
Jf Yy ? 0 then Y 4 is fathomed so go to Step la with 
PNT = j (Bound 10). If jJ = 1 and Y exceeds its upper 
bound from Bound 3 or NES 0, go to Step 10. Otherwise 
go to Step 5' with PNT = j. 
3) Step 5' With PNT = j determine capacity cost a = 
Cly,;). If a > € go to Step 4b! with PNT = N-l. Other- 
wise go to Step 6' with PNT = j. 
4) Step 6' With PNT = jJ, determine the related optimal 
production-inventory schedule by inspection and the re- 
lated cost PI(y;). Total cost is T(y 4) = Cys) + PI(yj). 
Go to Step 7' with PNT = j. 
ED peep (te With PNT = j if T(y4) > 23 go to Step 4b! 
with PNT = N-l. If Tlyj) = i. go to Step 9'. Other- 
wise go to Step 8. 
6) Step 9' Either include T(y4) - z on the list of 


current optimal solutions or replace the current list 





UM 
with m = T(y4), whichever is appropriate. Go to Step 
4b' with PNT = Nel. 


Monotone non-increasing demand 


Assume that demand is monotone non-increasing 
and that monotonicity is not destroyed by initial in- 
ventory, i.e. 

> 
d) =- Vo 7 d2 
By assumption (b) production should always take place 
as late as possible. Then for the CPI model the only 
optimal capacity expansion policy can be to expand 


capacity in the first period to satisfy the demand in 


that period, i.e. 
wW = dx - V -W 


This will then be sufficient to produce the remaining 
demands, so no capacity expansion problem need be 
solved. If the cost assumptions do not hold however, 
then this case is just handled with the general CPI 


model. 
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II.7 Transformation for P2 


A similar procedure to that used in trans- 
forming the CPI model Pl will be used to transform the 
CIDPI model Pe. First a lemma will be proved which 
will be useful in interpreting the transformed model, 
and in the algorithm to follow. 

Lemma 2.9 For any period i either 


0 


(1) w; > 0 and r, 


or 


l 
O 


(2) w4 Mand Ty = 
or 
(3) w = 0 and r} >O 
Proof: The proof is by contradiction and con- 
struction. Assume that Ws, > 0 and Lo O for any peri- 


od 1 gives a minimum cost. There are three Cases: 


Case li» We =r, > 0 


I i 
Let w! = w, - r, and r' = O 
T 1 4 
Then 
cCıw'-f,r'+t(w'=-r!)» N m, = C,(w,-r,)-f,°0+(w,-r,;)- N m 
D ie i: a 1 91 1 1 1 ^ e 


N 
= = + = e E 
CyWy- Cyryt(wy-ry) e 


N 
« C4W4- DILE LL] 


03. 


since na FE 


The assumption is contradicted and (1) has 
been constructed. 
Case 2 Ww. - r, = 0 

Let u = 0 and r! 


H 
O 


Then certainly 


Since cy > fy 

The assumption is again contradicted and (2) 
has been constructed. 
Case 3 w - ry < 0 

Let En - 0 and = r, =- Wy 


T 
Then 


N N 
cew'zt.p't(w;-r.)* E m, = CG,.*Ü0-f.(r4.-w4,)t(was-ra4)« E m 


N 
rien) Zum 


A 


N 
SUMUS ICI.) 5m) 


since Cy > Ps 
The contradiction is complete and (3) has 
been constructed. 


QED 
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Now using Standard Transformation 1 as in 
Section II.2, constraint sets (P2.1) and (P2.2) can be 


combined to give the set of relations 


UD SAYS = 2y] + Vo 


Ee) 
(wa-r3) = eet ys-ceyo * y] 


MR ea oN one) t Yy = 2YN-1 + XN-2 


Note here that, by Lemma 2.9 and Ww, r 2 0, if the right- 
hand side of a constraint in (2.11) is positive, then 

w > O and 1 0. If it is zero, then w, =r, = 0. 

If it is negative, then w, = 0 and Bi? 0. 


Solving (2.11) for fr and substituting in (P2) 


gives 
N 
P5: min Zg = an + hyvy + DE * (c4-£4)w.] + M 
S.t. 
X1 = SE yj] = "Vo 
(Pore) 


JA 


X4 dy + Y 4 - Maca 1 = ZEN 
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2,2) Xi t Vi l =V] KB. TESE os deese T md 


025,3) 
di * y, - y4.1 50 4 = 2,...,N 
Wy 2 dy * yj - (wg tvo) 
(P5.4) Wo = do - di * y5 = 2y) + vo 
ei 2? d4 - d, 4 + Yı - 2yı_ı + Yı-2 1 = 3,..,N 
x; Y, V, ¥ 2 0 and INTEGER 
where 
A AP A my ds 1,...,N-2 
ES = fn-ı - 2fy + My] - My 


D" 
ix TN " TN 


Nel 
M x TOME: 00 t Pydy + RAUCH = > Wo + (fo-fj-mi)v, 


Constraint sets (P5.1) and (P5.2) and the 
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first two sets of terms in the objective function are 
Just problem Pl! and thus are equivalent to a transporta- 
tion problem. Writing constraint set (P5.4) separately 
with the last set of terms of the objective function 


gives 


N 
P6: min Lg = ae wi 


> 
Det Y= Wo 


= 
t 


Iv 


Where ae Hay Eyy- + YN=22 


IV 


O and INTEGER 


which, given y, is solvable for w by inspection since 
the coefficients in the objective function are positive 
and the right-hand sides of the constraints are integer. 
Furthermore, it is possible that the right-hand side of 
some constraint i is negative or zero. Recalling Lemma 
2.9, it follows that in this case w = O andr, equals 
the absolute value of the right-hand side. So, for 
each possible y in the CIDPI model, it is necessary to 
consider a transportation problem and a trivial linear 


program. 





67. 

Essentially the same branch and bound enumer- 
ation scheme can be followed here as was used in the 
capacity increase case. Lemma 2.1 still applies as a 
finite bound on the enumeration and (2.5) is still a 
valid set of global upper bounds. 

The lower bounds on y are of a slightly dif- 
ferent form than (2.4) and can be gotten from (P5.3) and 


VES 0 as 


V 


Yous Yo. > C1 


> 
(2.12) Ya = Yı-ı - 44 ON 


11.8 Properties and Algorithm for P5 


The properties of P5 are similar to those of P3. 
Because of the new lower bounds (2.12) however, most of the re- 
sults must be restated and reproved. 
The statement and proof of Lemma 2.2 hold equiva- 
lently for P5 if the change of variables 
Wy = (w ~ r) ESLER NN 


is made. To arrive at sufficient conditions for convexity to 


hold in the optimal values of Z., the remaining results must be 


52 
relterated. For Lemmas 2.10, 2.11 and 2.12 and Theorems 2.5 and 
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and 2.6, when given a set of yi» IJ t is assumed 


that the remaining Y 54h are at their lower bounds for ke1,...,N-—j. 


Lemma 2.3 is restated as 

Lemma 2.10 If for given Yay 1 = ],...,J the lower 
bounds of Y 3+k are positive, k = 1,..., N-j then 
(W4 e k-T Jk) mO for k = 2,..., N=J. 


Prost: For any Ke{2,3,...,N-j} from (2.12) 


> 
Yj+k 7 Yj+k-1 ~ Ag+k > 0 


and 


> 
Yj+k-l ~ Yj+k-2 7 Qj4k-1 ? O 
where equality holds on the left. Then from (2.11) 


(wWj4ktj+k) = (YJok-YJ ke1*0 (Y J 4-1 Y 4 4k 22*0j 4 k-1) 
= 0 


QED 


Next,Lemma 2.4 becomes 
Lemma 2.11 If the conditions of Lemma 2.10 hold and 


y; = y,*n top mebbss o then 
(w,-r4)' = (W47r4) +n 


and for k = 1,2,...,N-j the lower bound on ux is 
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Yj+k tn 
and 
(44177443) ' = (W34+17P 547) -n 
and 
(W3+k=r3g+k)' = 0 A E 
Proof: From (2,11) 
AY = d; - di-1 + (y ,*n) = 2y 4-1 + Jj-2 
= (wi-r;) + n 

For k = 1 


> 


Es di+] = (y,*n) - d3+1 = Yj+l +n 


t 1 
Eu J 
so from (2.11) 

(44177 441)' = d3+1 = d; + (y 5+, +n) ~ 2(y,*n) t vie] 


= (wje Tje) n 


Doing induction on k, for k = 2 


' 2 ' = 
NON n ate 


(y441*n) 7 diio 





TO. 
Assume the relationship holds for a general k. 


Then for k+l 


t = t - d 
uc jae | tk 


(Yj+ktn) = djtk+l 


y + 
“j+k+1 
Finally by Lemma 2.10 


(8 eT 3 ec? ' =. k= 2,...,N 


QED 


The following is an additional property which will be 
needed. 
Lemma 2.12 If the conditions of Lemma 2.10 hold, then 
(w4,1-77441)' 7 0. 

Proof: From (2.11) 


1 1 1 < 
(n3 4177 444 7' A ($417 419344) = Ga - 0 


since Yl is at its lower bound so the first term is 
zero, and the second term must be non-negative by (2.12) 


QED 


Now Theorem 2.1 can be restated as 





TE: 
Theorem 2.5 If for given Y4> 1 = l,...,J, the lower 


bounds on Yj+x are positive for k = 1,...,N-j, then 


Cy, BS suc Peonvexein n £ermellj2,3,..). 


Proof: For y: - y,*n, ne{1,2,3,...}from Lemmas 


2.ll and 2.12 


(w,-7rj)' = (w 47r) dg 
(w 141774 43)' T (W441-7T441) - 1 s 0 


(WT quj)" = 0 KERSTIN 


The capacity change variables in period j could start 


from rs ^ 0 and go to Ws > 0 as n increases. Assume 
there is an ny such that E >=a0=2ror nue n4 and Wj > 0 


or n > nj. Then the change in capacity cost for each 


increment of n < ni 1s 
while for n 2 nj the cost change 1s 
SEUA A 


where A, > 0 follows from the cost assumption of P2 





12% 


that 


N N 


C, + > fj + Reik MJ. yel... N-1 


m PEE 
Further from P2 it was assumed that Cj > f} giving the 


convexity result 


QED 


Comet anvc2. 5.1. cita Is convex ldinin for metil, 2 
IA 
Theorem 2.2 and Corollary 2.2.1 still hold 
since, as has been pointed out, the production-inven- 
tory section of the model is equivalent. The main 
result, Theorem 2.3, must be restated as 
Bicone come lr eyhe conditions of Theorem 2.5 hold, 
then T(y , *n) ts convex in nakor nel 2 Bda 
Proof: The proof follows from Theorem 2.5 and 
Corollary 2.2.2 and the fact that the sum of convex 
functions is convex. 
QED 
Corollary 2.3.1 is unchanged. The above properties of 
the P5 variables and convexity in the optimal values of 
Z5 w111, with a few exceptions, permit the bounds and 


optimality tests of section II.4 to be used. 
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Algorithm for the CIDPI model 


The bounds and optimality tests of section 
II.4 hold for the CIDPI model with a few exceptions. 
Bound 3 is no longer valid since the possibility of 
decreasing capacity allows a Wy to exceed these bounds 
and not violate a. For Bound 4 the lower bounds are 
now found from (2.12). Optimality Test 3 can no longer 
use Bound 3. Finally, since Theorem 2.5 only gives the 
weaker result of convexity for C(y), Bound 8 has to be 
rewritten to ensure that the test is being conducted 
in the area where the slope of C(y) is positive. 
Bound 8' For given y,, 1=1,...,j if y, > 0,4 = eal, 
...,N, C(y;) = C(y,*1), and C(yj*1) 2 à then y; is 
fathomed for the preceding set of y4. 

The basic algorithm of section II.4 requires 
some changes for the CIDPI model. The lower bounds in 
Steps 3 and 4b are now found from (2.12) and the initial 


value for v in Step™3 is 
yj = max [(vg-d4), 0] 
The first part of Step 5 must be rewritten to reflect 


Bound 8' as follows: With PNT = j determine capacity 


cost a = CCy 4). If a > a and either j =N or j> N, 
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yi » u for i = J+l,...,N with C(y,) 2 C(y,-1) go to 
Step Ya with PNT = jJ (Bounds 8' and 9)..... The remain- 
der of Step 5 is the same. Finally, the update of 
Bound 3 in Step 9 can be dropped. 

The remarks of section II.5 pertaining to 
alternate optimal solutions are relevant here also. 
Finally, note that lower bounds (2.12) will serve to 
increase the feasible space for y and thus longer 


average running times should result for the CIDPI model. 


As with the CPI model, the efficiency of the 
algorithm for the CIDPI model is increased greatly when 
monotonicity exists in demand. The following hold for 
this section 

a) The cost coefficients of P2 (Chapter I) 

are monotone non-increasing with time. 


b) The cost assumptions of P2 hold. 
Monotone non-decreasing demand case 


Lemma 2.13 If assumptions (a) and (b) of this section 


hold, demand is monotone non-decreasing and 
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EL 
Nw, + Ver epi 
then there is an optimal solution with yy = 0. 
Proof: The proof follows trivially from the fact 
that these are the same conditions as for Theorem 2.4. 
Now, in addition to the possibility of eliminating 
idle capacity by the arguments of Theorem 2.4, there is 
the added possibility of scrapping it and recovering a 
salvage value, so the results hold a fortiori. QED 
The conclusions of section II.6 which follow 
from Theorem 2.4 are valid here for the result above; 
i.e. the transportation subproblem becomes trivial and 
can be solved by inspection for a given y, and a much 
more efficient set of upper bounds can be included 
implicitly in the algorithm. That is, Bound 10 of sec- 
tion II.6 can be used instead of Bounds 6, 7, 8, and 9. 
The algorithm changes of section II.6 remain valid with 
the exception that the new lower bounds in Step 4b' are 


found from (2.12). 
Monotone non-increasing demand 


For this case it must be assumed that 


> 


in order that the monotonicity be preserved. 
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Lemma 2.1H If demands are monotone non-increasing, 


assumptions (a) and (b) hold and further 
(c) f, + (Lom € cy MESE dE. (NZD 


then w = d4 - Wg - Vg and w, - DEM easy lis 
Proof: The statement for wj results from the fact 
that no backordering is allowed and by cost assumptions 
(a) and (b) and monotone demands, it can never pay to 
build extra capacity to produce early.  Assumption (c) 
eliminates the possibility of scrapping, capacity in any 
period to the extent that capacity must be rebuilt 
later, ,which, by monotonicity of demand, produces the 


result that w, = 0, 1 = 2,...,N. 


QED 


Lemma 2.15 If demands are monotone non-increasing and 
f 2 0 then MEE OF, 

Proof: The proof follows trivially from the fact 
that demand is monotone non-increasing so, with ape Us 
any idle capacity in a period can be salvaged for 
urcher revenues. 


QED 


As was pointed out in section II.6, when 





ES 
XI e 0, all capacity is used. The transportation sub- 
problem of P5 thus can be solved by inspection for a 
given y, speeding up the CIDPI algorithm greatly. The 
remaining remarks on algorithm changes for this case 
are the same as for the preceding case of monotone non- 


decreasing demands. 





Chapter III 


In this chapter, a special case of the CIDPI 
model will be presented which leads to a formulation of 
the entire model as a network. Other extensions of the 
basic model are then presented to allow piecewise- 
linear, convex production costs and discontinuities in 


capacity acquisition. 


III.1 Network formulation of the CIDPI model 

The formulations of Chapter II required a 
problem decomposition approach to provide suitable al- 
gorithms which took advantage of the problem structure. 
There is a special case of the CIDPI model, however, 
where the formulation can be simplified to that of a 
network and solved with a transshipment algorithm. 
This occurs when the cost of capacity increase is equal 
to the salvage value in any period. 

ouch a model might be useful when there is a 
large second-hand market in the units of capacity. A 
firm might then acquire or sell units with equal ease in 
any period at the going market price. Maintenance 
costs then reflect the investment needed to keep each 


unit at its rated capacity for possible future sale. 
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The derivation of this case comes from recon- 
sideration of the CPI model Pl. If WŽ 0 is dropped, 

allowing W to be unrestricted, the problem would still 

be bounded if capacity costs Cy, 1 = 1,...,N were mono- 
tone non-increasing, since from constraint (P1.1) and 


X 2 0 it must still be true that 


livre 


> a 
j WW = 0 i = JM SN 
This implies that each Ws is bounded below by the total 


capacity built up to period i, i.e. 
> 
Wi = m ZN; Mem y 


since if the costs are monotone, building and subsequent 
scrapping of capacity can never result ina profit. 

That is, with cy 2 0 and monotone i = IN, the 
right-hand side will always be as small in absolute 
value as possible. In fact, as a limit, it could never 
pay to build more capacity in any period than what is 
needed to produce all remaining demand in that period, 
since the excess capacity would be idle and can never 
reduce the total cost. Thus, with the assumption, which 


will hold for the remainder of this section, that 
> 


S 0 
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a set of lower bounds on W can be written as 


iv 


Wi “Wo 


S N 
Wo - Vo = cane! 
N 
Wd - L d 
i=jel 


tv 


i ml 3 ....,N 


Now with Y unrestricted in Pl, if Standard 
Transformation 1 is performed, then P3 no longer con- 
tains constraint set (P3.3). The integrality of y 
still holds, however. To transform the remainder of P3 
into a network formulation, slack variables 8 2 O must 
be added to constraint set (P3.1). Then the sum of set 
(P3.1) 1s subtracted from the sum of set (P3.2) to give 
a redundant constraint. Finally, multiplying the con- 
straints of set (P3.2) by -1, the transformed problem 


can be written as 
N 1 
P6: min Z¢ = T des tuer ciy4d* K 
Sets 


26.19 cy ype 5718 dy - v 


X4 + J1_] m J4 + S4 = dy i = 2 UN 





GIS 


-X4 + V -d4 + S 
(0% 2) AS di PA ESI 
-Xn - Yy-1 = -dy 
N 
(E63) 1 E D 


X, y, V, 520 and INTEGER 


where the constraint (P6.3) is redundant, the cy and K 
are as given in P3, and vy has been incorporated ex- 
plicitly in dy. 

The coefficient matrix associated with P6 
satisfies the sufficiency conditions for total uni- 
modularity as described by Hu (Hu 69). This guarantees 
that, since the right-hand sídes of the constraints are 
integer, the linear program associated with P6 without 
the redundant constraint (P6.3) will terminate integer. 
P6 posses more useful properties however, since the 
problem as written can be represented by the network 
displayed in Figure III-1. Thus the problem can be 
Solved using a transshipment algorithm to obtain the 
minimum cost flow. 

The only evident problems with the network 
are the cycles involving y and S which, if negative in 
cost, cause unboundedness. However, since flow S, is 


costless,the monotonicity assumption on the Ci is 
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S 








ee? 


Figure III-11 


Network for model P6 





Ithe figures on the source and sink arcs represent 


the arc capacities. For the remaining arcs the figures 


represent the flow variables associated with the arcs. 





eo 
sufficient to prevent this as seen from the following 
lemma. 
Lemma 3.1 If capacity costs = are monotone non-in- 


cnegsnpefore #7 .6.,N> then 


N 


gaye SNO JE TN SN 


Proof: It is obvious that 


> 
For j = Nel 


(ey-1-20N+My_73-My) + (eytmy) 7 (ey j-CN) * my. j Z o 


and by induction for j = k, l1 É k < N-l 


et slep cp) + m2 0 
QED 

It is obvious from Figure III-l that the 
Only cycles are of the form 

Yu» JN-1» ***» Jj» S ENS 1,95 N 
which have non-negative costs by Lemma 3.1. This con- 
firms the fact that, since P6 is obtained by a linear 
transformation on the bounded problem Pl, it too must 
be bounded. Emphasis must also be given to the fact 


that the capacity change costs 2r incorporated in P6, 


are completely non-stationary and thus admit to dis- 
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counting or any other aberrant structure. 


Networks and the general CPI and CIDPI models 


The existence of unimodularity in the special 
case of P6 permitted use of a network algorithm to ob- 
tain an optimal integer solution. The general models, 
Pl and P2, could also be formulated as networks but 
equivalence of the networks to the original models re- 
quires sets of constraints external to the networks. 
Moreover, these constraints violate the unimodularity 
conditions, as shown by the following examples, and 
thus exclude the use of any existing network algorithm 
to gain an integer solution. These examples are both 
linear programs having non-integer solutions with non- 
integer objective function values (eliminating alterna- 
tive optima). 

Example 1: Example for the CPI model 
Ror Neo moea K 20500519519518,18),-p.94(5,5,5,7,3; 
A Tee (3,3,353,2,2), di (33739, 
13,18,21), noe l= 1 
LP optimum 
EA AT GEA 7 A Ag” 15.25 
ve 6.25, vy = 8.5, Vo = 5.75 
A *.8.25 





00% 
O0OpJectIivye Function = 735,25 


Example 2: Example for the CIDPI model 
The problem is as in Example 1 with the exception 
of f = (8,8,7,7,6,6), d = (3,5,8,14,9,7) 
LP optimum 
A as ce Tb, xo» 7.0 
vı = 5.6, V2 = 8.2, v3 = 7.8, vy = 1.4, VUE 
w = 6.6 rg = .6 


Objective Function = 457.2 


III.2 Block Capacity 


In this section, CPI and CIDPI models will be 
extended to handle the case of block capacity. This 
terminology embraces any situation where capacity can 
only be acquired in units which are each capable of 
producing large numbers of units of demand in each 
period, as opposed to the implicit assumption in effect 
so far, that each unit of capacity produced a single 
demand unit per period. Block capacity, then, corres- 
ponds to the economists’ problem of "lumpiness" in 
acquisition which typically causes serious difficulty 
in a marginal analysis approach. 


In the remainder of this section it is assumed 





Do 
that in any period 
l Capacity Unit = k Demand Units 
where k is a positive integer and is stationary, and 
that initial capacity, Wo» is expressed in units of 


ELM 


CPI-Block capacity model 


Working first on the CPI model, in terms of 
block capacity, (2.1) ean be rewritten as 


> 
kw] -= Nes di we (kw +VQ) 


tv eee 


k (Nw +...+wWy) - W dy + ... + dy - (Nkw, *v,) 


N 


Following the same development as in Chapter II, and 


using Standard Transformation 1 to solve for w gives 


[Yi + di ze (kwotv.)] 


[Jo - 2y4 * do - di * vg] 


z 
ho 
y 
SP ir =/1r 


[yg 7 ?3y.3 * Xy. 2 * dy 7 dy. 





^Depending on the source of the demand, management 
has the option of choosing the most appropriate unit for 


demand, e.g. box, pallet-load, car-load, etc. The 
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and the expressions for w, can now be substituted back 
into Pl to get the transformed problem. In this case, 
substitution in the W 2 0 constraints gives somewhat 
expanded results. The constraint set corresponding to 


(P3.3) now looks like 


T > 
[yy + di ~ (kwotv,) J O and INTEGER 


- [y, - 2yg., * dy - dy ,] ? 0 and INTEGER 


Taking the first constraint and ignoring the integer 
requirement gives 


> 
Mi KW + View = di 


while the integer requirement gives 


l 
s Ly, 4 di =- (kwotvo)] - n4 


where n4 is a non-negative integer. But this is the 


Same as 


= - + + 
y nik d KW V 





value of k can then be set appropriately. 
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A similar development for period j gives 


2 
SA ya d^ 


and 


ae metta s 


for n, 2 0 and integer. So, in general, substitution 
in W 2 0 gives two sets of constraints, the first sim- 


1lar to set (P3.3) 


tv 


l O O 1 
(3.1) a> +d dn =V 
: Yo = “Vy ae eG 
> 
y 20 


and the second, a set of integrality constraints on y 


Yi = nk - d4 t KW + de 
i - - z 
Rue er ots n, 2 0 and INTEGER 
AN 
yy" nk + Yx-1 - Yy_o + dy 1 = dy 


The effect of k, then is to superimpose a lattice of 
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interval k on the space of y. Up until now, k = 1 and 
so the enumeration on y had to search every feasible 
integer point in y space. As k becomes larger, however, 
the number of enumerable points in y space becomes 
smaller for a given d and the problem can be solved 
more quickly. 

Doing the remaining substitution for w in Pl 
gives the following form of the CPI-Block capacity 


model. 


n mn, E i [p,x,thsv, ] + zi 3 c'y +K] =- (k-1)6 w 
To a1 134 1157 paeii me 


Sauce 


Coors Grsee ys (sel), and (3.2) 


0 and INTEGER 


iv 


X, V 


<i 
Iv 
O 


where ER and K are as in P3. The basic algorithm is 
valid for P7 with a few changes. In step 4b let 
Y Yj + k. The global upper bounds are changed 
slightly by the parameter k, since the maximum value 
for w] is now 

1 


N 
W, = =-( IE d,-v,) -w 
l mc e O 
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So using (3.2), (2.5) of Chapter II becomes 


3.3) yj S jkemax Wi - t d, * v 


dT Yo 0 d t Lan 


Nngorder eto satisiy (3.1), (3.2) and y 2 0, the lower 
bounds on y are first obtained from (3.1). If y 2 0, 
then (3.2) is also satisfied. If y, « 0,15 1$ N, 
then by adding m, K to Ys for successive values of 


ee a SLI 2 0, a feasible lower bound 


can be obtained which satisfies (3.2). 
CIDPI-Block capacity model 


Now for the CIDPI model, (2.11) can be re- 
written in terms of block capacity as 
n 
(w1-r1) a Lay + Ji mE (kwo*v )] 
i 
(w5-r5) = 855 - d) ~- 2y; +t yo tł Vl 


H 
1 = outs 


Substituting for r 2 0 in P2 glves 


d » 
Wi - zm + yı - (kwotv,) J - 0 and INTEGER 
« 
4 





ze 


2 > 
Wy - LAN - dy-1 + Yy” 2YN-1 + yN-2J - 0 and INTEGER 


Following a development similar to the first part of 
this section, the zen constraint is equivalent to the 
two constraints 
l 
> 
WO eyed jetty cey3 1*3 3.2) 
and 


Y 4 = nk + dj} = dj + 2Y 3-1] = J j=? 


where n; is any integer, but now a < 0 is allowable 
since the right-hand side of the first constraint is 
no longer required to be a non-negative. 

In general, the two constraint sets which 


result are 


tv 


i 
Wy LA ty yr Wee): 


tv 


1l 
(3.4) Wo rLda-dj-ey*ystv, J 


Vv 


jJ 


and 


M nik - d4 + (kwotvo) 
M) 
Yo = nok + ey] * dj - d5 - Vo 
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Yi + dn- d n 


o 3 INTEGER 


= r 
Y3 n3k ey» 1 
O UN 


Yu ? nyk + 2Yy-1 - Yn-2 * dy] - Ay 


Substitution for r in the remainder of P2 Just repeats 


the first three constraint sets of P5 and gives the 


model 
P8 in Z 2 [ thav,+f" ( ff) wees M' 
: min = pax V y t(c,- W + 
A iA 
Sata 
CES. len oe) cand), and (3.5) 
*, V, W 20 and INTEGER 
y20 
where 
f" = ip 
i k i 
M! s iu - UD y 
k k l o 


Dor m and M from P5. 

Modifications to the basic algorithm for P8 
are the same as for P7 except that, if the lower bound 
for any y, is positive, then (n,;k) should be subtracted 


from y, for successive values n, = 1,2,3,..., until the 
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smallest y 2 0 is obtained. This eorresponds to the 


unrestricted n in constraint set (3.5) which is caused 
by the possibility of having Rc Urn any perdos 
It should be noted here that the special 
cases of Chapter Il are valid here only for k = 1. 
Another point worth noting is the increase in danger of 
rounding off from the related linear programming solu- 
tion as k increases. Not only are the unit capacity 
costs larger for large k, but also shifting a unit from 
one period to another has a proportionally greater 
effect on the optimal production-inventory schedule. 
Finally, the increase in the lattice interval 
with k leads to a useful scheme when deciding on a capac- 
ity increase (and decrease) strategy. If several sizes 
of capacity are available, then the algorithm can be run 
for each value of k with the running time probably de- 
creasing for each larger value of k. This allows a third 
dimension of optimization for management with the cost 
of gaining the added information governed by the algo- 
rithm running time for the smallest value of k. Note, 
however, that the form of Standerd Transformation 1 does 


not permit k to be non-stationary within a model. 
III.3 Convex Production Cost 


Up until now, it has been assumed that all 
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production was being carried out by one shift working 
regular hours with unit production costs p. In this 
section the CPI and CIDPI models with block capacity 
will be extended to allow multi-shift and overtime pro- 
duction. Although the work done here is for one addi- 
tional shift, further extensions are trivial, as will 
be pointed out later. 

The crucial assumption of this section is that 
the unit cost of production increases for each shift 
and for overtime. In the manufacturing context, this 
assumption is not restrictive and actually models real- 
ity, since it is commonly accepted that there exist 
"shift differentials" in manufacturing costs for a 
second and third shift. These can be ascribed, in part, 
to the increased time spent on maintenance and to 
greater human error in succeeding shifts. As for over- 
time production, the wage differential alone is suffi- 
cient to justify tne assumption. 

Let 
ie second shift production in period 1 
s second shift unit production cost in 
period i 
Production capacity has been defined in terms of regu- 
lar time production. It is assumed that, in terms of 


the second shift, an equal amount of capacity is 
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available. Thus 


< 
? = - 
iM Te + Wo i NN 
Assume that 
Py = P} i = N 


Then the production cost function for period i 


' = tyt 
an) TES 4 DD 


is plecewise-linear convex. Similarly, with a suitable 


rearrangement of terms, total production cost 


N 
P= E 


Pex. t) 
j=] $ i? 


is piecewise-linear convex. 


CPI-Block capacity- Convex production cost model 





First the extension of the CPI model with 
block capacity will be presented. The term Te 
must be added to the objective function of P,, while 
the following constraints are added to the constraint 


set: 





202 
x' = k I w S kw A A 


Gee la) E 
27220 and INTEGER 


and the appropriate block capacity coefficients are in- 
serted in the remainder of the model. Then (P1.1), 
(Pl.1a), and (P1.2) can be combined to give a set of 
necessary and sufficient conditions for feasibility 


which take into account two-shift production. 


> 
ekw, S di = Yo = ekwa 


> 
2k(2w two) VE ss di + d; - Vs 7 li kw 


N 
2 


7 Nee fet o 


Again by Standard Transformation 1 a solution for w can 
be obtained. Then substituting into w 2 0 of Pl gives 


the following two sets of constraints 


> 
Y4 x t ekw - di 


iv 


(3.6) Yo 2x Tg. = d = Vo 


L 2 


> = 


and 
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yi z en,k - d, + 2kWo + ME 
y» == enok + ey] + di = d, - Vo 
p 2n-k + 2 d d = 
y = n t y - y + - ny - 
3 3 € 1 e 3 and INTEGER 
e La 


INTE enyK * 2yN-] - YN-2 + dyiy - dy 


From (3.7) it is clear that the lattice interval on the 
y space has been increased to 2k with the introduction 
of a second shift. 

Completing the substitution and transformation 
from Pl gives 


N 
P9: min Z, = à X4 * p'x' + hv. + c*y ] + k! 
Tc genet PX, et £4 


Set. 
< 1 
2 = 5 di + Y y Vo) 
xa - ss + y, P y4 3) 1 = c,...,N 
yos + = 

x (d, yi 7 V ) 
(P9.2) A « 1 > 

x gd, + yy cv P diii 
(P9.3) x, + x! tvi- —- VQ *. d, 10221 ur," 


(3.6) 
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(3.7) 


2 0 and INTEGER 


- 

wv 
> 

vw 
«I 
i 


tv 


y 0 
where 
" oz E t zm 
and 


K'= 5 K - GE. we 
for e. and K from P3. 

The addition of the extra term to the ob- 
jective function and the extra set of constraints 
modifies the transportation problem structure of the 
subproblem. In effect, N additional sources have been 
added and the number of arces has been doubled. This 
development might make it beneficial to reformulate 
the sub-problem as a transshipment problem which can 
be solved with an appropriate minimum cost flow algor- 
ithm (Ford-Fulkerson 62). An investigation would have 
to be made of the tradeoff between the faster algorithm 
but higher network growth rate of the transportation 


problem against the slower algorithm but lower network 


growth rate of the transshipment problem. Other algo- 


rithm changes and remarks on the bounds for y are 


2 
identical to those for the CPI model of the previous 


section. 


CIDPI-Block capacity- Convex production cost model 


The convex production cost case of the CIDPI 
model with block capacity can be obtained easily by a 
similar development. In P2, the term 3 H is added 


i21 
to the objective function, the constraints 


i 
(P2.la) x' - k E (w;-r;) = kw al oe eee 


O 


Y! 


X' 20 and INTEGER 


are added to the constraint set and the appropriate 
block capacity coefficients are inserted in the remain- 
der of the model. A set of necessary and sufficient 
conditions for feasibility is then obtained from (P2.1), 


(P2.la), and (P2.2) as 


> 
2k(w4-ri) SRM- do- vs 2kw, 


l 
ek[2(w)-r])+t(w,-r,)] N d, Dados VR li kw e 
CICER a a 2Nkw 


Using Standard Transformation 1 to solve for r and 
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substituting for r 2 0 in P2, the following two sets 


of constraints, analogous to (3.4) and (3.5), are ob- 


tained. 
il: 
> 
Wy - Std tyr -vocekw,] 
(3.8 Z ENS d,-2y.*v ] 
3:0) Wo = 5 9794-29, s 
a 
BE une) OO 
and 
y, = en,k - di * E. + ekwa 
Y. * enok + ey, + d4- da = v 
(3.9) 2 > = 


= r = + = n, INTEGER 
en3k ey, Y, d, d3 1 


d BEIM o ePi 


Yy * 2Myk + 2YN-1 - Yyz2 * dg. - dw 


Again note that the lattice interval on the y space 
has been increased to 2k. 
Completing the substitution for r in P2, the 
model is transformed to the following: 
E = N 1 1 
P10: min 244 deroga tp seh 


n n 
D a a +M 





LONS 


Sie bs 
l 
« 
Xy - Mil + Yi - ME 
(P10.1) l 
x, Í Í(d, + x ) mer N 
re 211 u: 
(£1 " 
x 5 (dy UM v ) 
> a BU + Ys — Y4-}) 1 = P NEAN 
(P10.3) T T x. t Ne | TX Va = d, I = B. 
(P10.4) 5 
elie © LS DN 


(3.8) 


(3.9) 


2 0 and INTEGER 


ot 
w 
ex d 
he 
<i 
e 
= 4 
i 


«Ii 
iv 


where 


1 
"n - = 
n I i Pe EMEN 
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and 


where ir and M are from P5. 

The remarks of the previous section relevant 
to the changes in the algorithm and in the bounds apply 
here also. As noted there, the subproblem can be form- 
ulated as either a transportation problem or a trans- 
shipment problem, 

As has been mentioned, the techniques of this 
section easily extend the models to encompass a third 
shift and overtime work. The capacity available for 
overtime can be assumed to be equal to that for one 
regular time shift. This presents a problem since 
overtime usually includes both time-and-one-half and 
double-time wage periods. Some subjective unit pro- 
Guction cost would have to be chosen to represent the 
entire period rather than differentiating between the 
two cost levels. 

Finally, as was already remarked, an add- 
ition of another shift increased the lattice interval 
on the y space. If the extension is continued, then 
the first term on the right-hand sides of (3.6) and 
(3.8) becomes 3nk with a third shift, and link with 


overtime, increasing the interval by nk each time. 





Chapter Ty 


The basic idea of the enumeration algorithms 
presented is similar to work done by Murty (Murty 68) 
and Gray (Gray 71). A discussion of these related 
approaches opens this chapter. Also included are a 
short discussion of the restart capability of the 
algorithms and a presentation of computational results. 
Finally, some conclusions regarding the preceding work 


and suggestions for possible extensions are given. 


IV.1 The Decomposition Approach to Integer Programming 


Since the algorithms developed in Chapter II 
use the idea of decomposition, it would be useful at 
this point to put them in perspective relative to some 
similar approaches in integer programming. This is not 
done to claim extensions to existing theory, but rather 
to show contrasts and similarities occuring when a 
common approach is used on different problem structures. 

One approach which uses the properties of 
duality theory to develop a gradient algorithm for 
mixed integer programming was given by Benders (Hu 69). 
A drawback to the algorithm is that efficient use of the 


gradient information requires the solution of a series 
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of integer programs where the number of constraints is 
increasing by one at each iteration. The algorithms 
of Chapter II only make partial use of the available 
gradient information in the convexity properties which 
are developed. Possible further use of this information 
was not investigated in the course of this work. 

Other approaches more closely related in 
principle and in spirit to the algorithms of this 
thesis are the work of Murty (Murty 68) and Gray (Gray 
71). Both works deal with an exact solution to the 
fixed charge problem and, as pointed out by Gray, the 
two have mutually enhancing properties and could pro- 
bably be combined to give a powerful algorithm. 

Murty develops a decomposition approach which 
ranks the extreme points of the convex hull according 
to the value of the variable cost. An upper bound on 
variable cost is given by the use of a minimum value for 
the fixed cost ina manner similar to that which gives 
an upper bound on capacity cost in the algorithms of 
this thesis. The variable cost upper bound can then 
be used to terminate the enumeration of extreme points. 
The algorithm's efficiency is greatest when the fixed 
costs are small both in absolute and relative magnitude 
with respect to the variable costs (Gray 71). 


Up until now, no mention has been made of the 
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fixed charge aspects of Pl and P2. It is not hard to 


make such an interpretation of the capacity costs how- 
ever, since once incurred, they are spread over the 
units produced in all future periods. The capacity 
cost of this thesis can thus be thought of in terms of 
a fixed charge although the problem formulation does 
not make this explicit. 

In light of the above, the approach of this 
thesis is much closer to that of Gray than of Murty, 
since Gray's decomposition approach works mainly on 
tne fixed cost section of the problem and deals with 
the variable cost as a subproblem. His fixed cost 
section is concerned with the optimal profile of 
open and closed routes of a transportation problem 
where each open route incurs a fixed charge. The scheme 
used is enumeration in a number system of base 2N (where 
n is the number of sinks) which takes advantage of both 
the transportation problem structure and the 0-1 prop- 
erty of the integer variables. In Gray's approach, an 
upper bound is developed on the fixed cost which is 
analogous to a of Chapter II. 

The major difference between Gray's (and 
Murty'S)model and those of this thesis is the 0-1 
nature of the integer variables. Transformation of Pl 


and P2 to 0-1 variables could be done, but at the 
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expense of a huge increase in the number of variables. 
Gray also encounters another difficulty not 
present in this thesis. He has no means of ensuring 
that each node of the enumeration is feasible to the 
to the main problem, thus expending resources in enumer- 
ating, and testing for, infeasible nodes. As has been 
mentioned previously, enumeration on y assures feasibil- 


ity at every node in the algorithms of this thesis. 


IV.2 Restart Capability of the Algorithm 


The algorithms developed have a feature that 
enhances their flexibility and usefulness, that is the 
capability of starting them at any point in the enumer- 
ation. This feature is inherent in the method of 
enumerating on y; since given any y' and the value of 
the algorithm pointer (PNT) plus some value for the 
current bound on capacity cost (a) and the current op- 
timal cost, the algorithm can be started from Step 5 
and run normally. Even if a and current optimal cost 
are set to infinity, the next execution of Steps 6 
through 9 will establish a new set of. bounds, though if 
the optimal solution precedes the point at which the 
algorithm is restarted, the new bounds may not be as 


tight as those which would have resulted from the 
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optimal solution. If the optimal solution is given by 
some vector after y it will be found whether or not 
initial bounds are provided, 

This restart capability, when combined with 
the linear programming bound, provides a very useful 
tool to management, especially when a large problem 
which could be expected to have a long running time is 
encountered. Instead of starting the algorithm from the 
lower bounds of y, a y' corresponding to a w slightly 
lower than any rounded-off solution to the linear pro- 
gram can be chosen as a starting point. The algorithm 
can then be run for a sufficient period of time to in- 
spect all w which could have been rounded-off from the 
linear programming optimum. An optimal solution is 
assured here if it equals [ts . If none such is 
found, then a management decision can be made based on 
the difference between the current optimum and Z* as 
to whether the full algorithm should be run in 2 
attempt to achieve further improvement. 

It might also be advantageous, e.g. because 
of limited computer time, to run the algorithm for 
relatively short periods of time. If no terminal 
condition has been reached then, again, a management 
decision can be made on whether to attempt to improve 


the current optimal solution with respect to a 





108. 
If so, the algorithm can be started right where 1t last 


terminated with no loss of efficiency. 


IV.3 Employment Smoothing 


The relevance of the literature on employment 
smoothing was mentioned in Chapter I. Now that the de- 
velopment has been completed in the context of capacity 
expansion, a reinterpretation of the results can be made 
to handle the employment smoothing problem (O'Malley 
et.al. 66, Silver 67, Lippman et.al 67a and 67b, Lippman- 
Yuan 69). Reference, in this section, will be to the 
most general models developed, i.e. those of section 
III.3, from which all special cases can follow trivially. 

The time period for the employment smoothing 
problem is much shorter than that for capacity expan- 
sion by the simple fact that hiring, firing, and usually 
even training, labor does not take as long as changing 
fixed capacity. It 1s assumed, therefore, that fixed 
capital capacity is never a constraint on the size of 
the work force. 

The work force units can be either in terms 
of individuals, or groups such as gangs, crews, etc. 

It is assumed that there is a capacity associated with 


each labor unit which can be expressed in terms of the 
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consumptive factor bundles 9,. This capacity could be 
gotten from the associated fixed capital optimal oper- 
ating level capacity or by a separate marginal pro- 
ductivity level argument for labor as was done for 
capital in Chapter I. It is again assumed that the 
production function exhibits constant returns to 
scale in labor and the consumptive factor. The accon- 
panying assumption is that the existing labor force can 
be used at any level up to its capacity with no lost 
opportunity cost of idle labor other than wages. 

A change in the factors of capacity will have 
to be made. Before, the capacity definition was ex- 
pressed in terms of one shift doing regular time work. 
Now capacity is defined in terms of all possible regular 
time output which can be achieved by the existing work 
force. Thus it is assumed that the shifts are filled 
sequentially, i.e. before adding second shift labor, a 
full first shift must be operating. 

The CPI model P9 would have limited appli- 
cation in the employment smoothing sense because it 
ignores capacity reduction. Its simplifying features 
relative to the CIDPI model however, make its use at- 
tractive when the circumstances permit. Two instances 
when such use would be feasible are those of a firm 


phasing in a new production line where the phase-in 
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period is relatively long, say a year; and that of a 
firm, with a reduced work force on an existing line, 
facing a projected growth in demand, 

Reinterpretation of P9 to employment smoothing 
leaves the inventory variables V and their coefficients 
h unchanged. The production variables X now represent 
all regular time production, and no differentiation in 
production costs p among shifts is allowed. Overtime 
production is x' and again undifferentiated overtime 
production costs p' must be assigned. The capacity ex- 
pansion variables w now become work force hiring var- 
jables and the associated costs c can include all the 
normal costs associated with hiring personnel, i.e. ad- 
vertising, personnel staff costs, training, etc. Wages 
now take the place of maintenance costs for m. 

As the reader has probably noticed, an inter- 
pretational difficulty exists concerning the capacity 
ascribed to overtime production in P9. The assumption 
for that model was that overtime capacity was approx- 
imately equal to the capacity of a regular time shift. 
What is necessary now, if Standard Transformation 1 is 
to be used, is the seemingly ludicrous assumption that 
overtime capacity is equal to all regular time capacity. 
With a time span of one to one and one-half years, how- 


ever, discounting will not be significant so that over- 
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time production costs will be much higher than any reg- 
ular time costs over the period. Furthermore, hiring 
costs are now in the same order of magnitude as pro- 
duction costs so that, within a reasonable range, 
hiring additional regular time labor should be cheaper 
than working any significant amount of overtime. A 
reasonable approach to the problem would be to first 
solve it for regular time only, and then try a solution 
permitting overtime. If the overtime assumption causes 
difficulty, then either a subjective modification of 
costs can be made and the problem re-solved or the 
approach has to be abandoned. For the majority of 
cases it would seem that the approach is useful. 

For the more realistic CIDPI model P10, the 
only additional interpretation necessary is that the r 
now represent the firing variables and the f the firing 
costs. The same difficulties exist for overtime capaci- 
ty and the previous remarks apply equivalently here. 

An argument can be made for the equivalence 
of hiring and firing costs in any period. Certainly 
what is included a..d excluded from each cost is sub- 
jective. Assuming that equality exists and, further, 
that the ¢ are monotone non-increasing with time, then 
the special case of the network formulation is appli- 


cable to the employment smoothing problem with the 
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preceding interpretations on the relevant variables 
where in addition Wy Oe corresponas, to the firing 
variable ae O of P10. Moreover, use of the network 
excludes the possibility of overtime and requires that 
k = 1, No obvious extension of the network model to 
remove these restrictions exists. 

Finally, the assignment of a value for k in 
P9 and P10 will undoubtedly be more difficult in terms 
of labor than it was in terms of fixed capital. On the 
other hand, the sensitivity of the models to k should 
not be significant for labor because the hiring and 
firing costs are much smaller than the corresponding 
capacity change costs, and thus the impact of fluctu- 
ations in the labor force level will not be as great. 

This work, then, becomes the first exact 
integer approach for the employment smoothing problem.? 
The algorithms given for the special demand cases in 
sections I1.6 and 11.9 are more general than those of 
Lippman et.al. (Lippman et.al. 67b), and algorithms to 


handle general demands and general linear costs have 





las was mentioned in Chapter I, a "round-off" 
approach from a linear programming model has been pre- 


sented (O'Malley et. al. 66). 
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been given in addition. Finally, the Similarity of 


the two areas of capacity expansion and employment 


smoothing certainly suggests that a combined effort 


would be useful. 


IV.4 Computational Results 


The algorithms developed in this thesis were 
coded in Fortran IV and run on the IBM 7094 in the 
Johns Hopkins University computing center. The test 
problems used were set up by the author using non-random 
data. Twelve problems were generated in all, each one 
as a ten-period problem. When run, however, each prob- 
lem was run eight times using the first k periods of 
data each time for kKi*# 3,...,10. Thus e.g. for problem 
10, the three period problem is numbered 10.1,while the 
ten period problem is numbered 10.8. This numbering 
system is followed in the tables of data which follow. 

The following conditions are satisfied by all 
of the problems: a) Costs are monotone non-increasing 
with time. In addition, the costs of Problems 7, 98, 9 
and 10 are stationary. b) k= 1 (l Capacity unit = 1 
demand unit) c) Only one shift of regular time work is 
being used. 


A linear program was solved for the CPI and 
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CIDPI case of each problem (using models Pl and P2). 
Because of the "nearly" unimodular structure of the 
general models, a large number of integer extreme points 
probably exist in each convex hull. In fact, for the 
166 linear programs solved for Problems 1 and 3 to 12, 
a total of 40 terminated integer. 

It had been thought that the number of periods 
would be the most significant factor in problem running 
time and that the appropriate curves could be presented. 
Instead, the models turn out to be sensitive to both the 
cost and demand structures of the individual problems, 
as well as being dependent on the number of periods. 

The computational results are, therefore, presented in 
tabular form in Tables IV-1, IV-2, and IV-3 at the end 
of this section. 

An example of the sensitivity of the algorithms 
to demand and cost structure is given by Problems 9 and 
l2 using the CPI algorithm. The data for the problems 
is as follows: 


Problem F 


Ci 


EMG ea rhe A O 13,195 16) 
( 30, 30, ..., 30) 

(25035. 0, 23) 

eas «2-5 5) 

au el) 


OI B CI 
| u I 
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Problem r2 


HERE EUSEB T. 9. 12, 10) 

c = (250, 240, 230, 220, 210, 200, 190, 
200,170, 160) 

ms (SIMS... 8) 

PI A O, 3) 3) 

DICE NES 2) 

For t = running time 
n = number of periods 


a curve fitting program was run for a polynomial 


an? 


et 
il 


and an exponential 
bn 


fit on the data. For Problem 9, using 7 data points the 


results were: 


Polynomial 


a = .00028 b = 5.10 

Coefficient of correlation = .960 
Exponential 

a = ,00635 b= .94 

Coefficient of correlation = .606 





qub. 
while using 7 data points for Problem 12 the results 


were; 


Polynomial 


a = .00006 b = 5.63 

Coefficient of correlation = .875 
Exponential 

a = .00573 b = ,87 

Coefficient of correlation = .973 


When the data for 6 common points was averaged the 


following resulted: 


Polynomial 


a = .00003 b = 6.21 

Coefficient of correlation = .561 
Exponential 

a = .00339 b= .993 

Coefficient of correlation = .770 


Similar results were obtained using the data from the 
other problems for both the CPI and CIDPI cases. 


The only conclusion which can be drawn, 


then, 


is that the number of periods involved 


nificant factor in problem running time, but 


exist other significant factors. 


Inspection 


tables will verify the claim made at the end 


II.8 that longer running times should result 


CIDPI algorithm. 


LIF. 
is a sig- 
there 
of the 
of section 


for the 


Two interesting examples remain to be shown, 


the first demonstrating the acknowledged fallacy of 


rounding-off the linear programming solution, and the 


second a demonstration of the power.of the restart ca- 


pability discussed in section IV.2. 


Example 1: 


able 


gramming optima respectively are 


For the CPI case of Problem 12.2, the vari- 


w 
^y 
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(with w 


3.25 
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V3 = 
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vy 


values resulting from the linear and integer pro- 


l and van = 1) 


O 


ad Oe Vy= 0 


The difference between building all capacity in period 


l and building it all in period 2 could be significant 
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to management. Furthermore, the production-inventory 
policies are quite different in each case. 
Example 2: The CPI algorithm for Problem 8.8 took 
27.067 seconds to terminate where termination resulted 
from the linear programming bound. The optimal capacity 
expansion vector for the linear programming solution 
was 


Bea 3815, DO «se, 0) 


The CPI algorithm was started with a y corresponding to 
Wi = 1, Wo = 4, ey 3, and the remaining yy through Yio 
set at their lower bounds. The time to termination by 
the linear stesse nme bound from this starting point 
was 5.9167 seconds. 

In terms of core usage in the computer, the 
algorithms are limited only by the requirements of the 
transportation problem. The version of the transporta- 
tion algorithm used requires only one mxn array where 
for the algorithms of ChapterII,m = n = N, while if 
a transportation algorithm is used for the convex pro- 
duction cost use of section III.3, m = ZN, and n * N. 
The enumeration itself is performed using a single 
vector N+2 words in length. The growth in core require- 
ments is therefore dominated entirely by the growth of 


the transportation array. 
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Table IV-1 


Computational results. for the CPI algoritnhm 
with monotone non-decreasing demand 


Periods | Optimal Values Run Time(sec) [Step (€ 
Number (N E N ! P Termination] 2 
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| 586 3.50 85.233 | No 
942 SIE 790350 No 


enumeration was stopped. 
“Step 8 uses [zr plas a termination rule for the 


enumeration. 


Table IV-2 


Computational results for the CPI algorithm 
with general demand 


roblem| Periods| Optimal Values | Run Time(sec) [Step 8 
E ac leu 
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9.1 
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2005339 exon O 
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Table IV-3 


Computational results for the CIDPI algorithm 
with general demand 


S| Optimal Values | Run Time(sec) [Step 8 


e Te el 2 unen 


302 57 1.067 | 
387 : 033 No 
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BeOS. 
ive Conclusions 


The algorithms developed in this thesis take 
advantage of the structure inherent in the combined 
capacity expansion and production-inventory problems 
when no backordering is allowed. The decomposition 
approach of handling the fixed-charge (capacity ex- 
pansion) and variable cost (production-inventory) por- 
tions of the problem separately allowed the application 
of special techniques for each part which could then 
be combined into an efficient algorithm. The regular 
behavior of the y variables in the fixed cost por- 
tion of the models allowed the development of useful 
optimal cost convexity properties in the variable cost 
section and in the overall problem. 

The main problem of integer programming al- 
gorithms has not, unfortunately, been avoided in this 
work. That is, these algorithms display the same 
unpredictability in running time displayed by all 
known integer algorithms., Furthermore, computational 
experience has shown that the algorithm running times 
are sensitive to demand and cost structures as well as 
being dependent on the number of periods involved. 

The flexibility provided by the restart 


capability of the algorithms is one of their most 
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attractive features from a management point of view. 
The feature allows for efficient use of the linear 
programming solution initially, by permitting a 

search among all rounded-off values. Available com- 
puter time can then be used efficiently for the main 
algorithm by running the program for short periods of 
time, since no information is lost when the enumeration 
is stopped. 

The economic problem of lumpiness in capacity 
acquisition has been handled rather easily in the 
emsdonsectschaptepn lil. In- -fact the more severe 
the lumpiness is, i.e. the higher the value of k, the 
larger the intervals in y-space where the enumeration 
is performed. 

Finally, the observation of an equivalence 
between the CPI and CIDPI models and the employment 
smoothing literature is of no small significance. The 
amount of work done in employment smoothing has been 
fairly extensive and unification of the results of 
these two areas would enhance the existing results in 


Boch» 


IV.6 Extensions 


The restrictive assumptions which led to the 





les, 


results of this thesis provide for some obvious ex- 
tensions. Perhaps the one with the greatest usefulness 
in the management context would be a provision for 
backordering, where the backordering horizon could be 
necrmuctedaco M perlods, for Me{l,2.3,...}. Unfortu- 
nately, incorporation of backordering in the model will 
destroy the structure of Standard Transformation 1 and, 
in addition, if a decomposition approach is taken to 
solve the problem, the size of the enumeration space 
would be greatly enlarged. On the surface it seems as 
if an entirely new approach must be developed for this 
extension., 

Another useful extension would be the develop- 
ment of a model which could handle concave costs of 
capacity expansion, i.e. decreasing marginal cost in 
che adailtion Oofæcapacity in any period. The optimal 
policy results of Manne and Veinott for capacity ex- 
pansion with concave costs (Manne-Veinott 67) might 
bezuserulshere, 

As was mentioned in Chapter I, the work in 
capital budgeting is just one step up in aggregation 
from the content of this thesis. A simple extension 
in the context of the work already done would be to 
have a capacity expansion budget in each period. This 


adds N budget constraints to the models and these 
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ponmstprasuUsceoudbldobe-chandled easily in the fixed charge 
(capacity) part of the enumeration algorithms. 

An extension to many production lines with an 
overall capacity budget in each period would not 
be so easy, even if each line had a single product. 

The difficulty will spring from the problem of opti- 
mally allocating the capacity resources among lines in 
a period as well as optimizing for each line over N 
periods., A very large-scale integer program will cer- 
rary result. 

Amma ler Scale put equally interesting. 
problem would be an extension to multiple products on 
Mice sane production Line. As was noted in Chapter TII, 
the work of this thesis limits k to a single value for 
any model because of the form of Standard Transformation 
1. Multiple products are sure to imply different pro- 
duction rates and, hence, different values of k. The 


challenge appears interesting. 
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